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ON SEISMIC COMPUTATIONS, WITH APPLICATIONS, IT 


M. M. SLOTNICK? 


ABSTRACT 

This paper is a continuation of a previous one given, with a somewhat similar 
title, in which are discussed the travel-times and other similar properties of the wave 
paths in a medium in which the seismic velocity is an integrable function of the depth. 
In this paper we consider similar properties of the waves from a point on the surface 
to a point in a “well” in such a medium. 

In this paper, it proposed merely to continue the results of a 
previous paper*® on the subject of wave paths and time-distance 
curves for media in which the seismic velocity is an integrable func- 
tion of depth. The “time-depth curves’’; i.e., the curves indicating 
the travel-times of seismic waves from a point of the surface to the 
points in a “well” are obtained in terms of ‘‘dimensionless” units in 
the so-called “‘linear” and “exponential’’ cases. It is hoped that these 
two cases, explained in some detail in this paper and in that cited 
above, will serve as a sufficient outline for other cases that may be 
proposed. 


THE LINEAR CASE: 9=2+ah. 


In this case, the wave-paths are the circles of the family whose 
equation is: 


(1) {a — [(z — pmo?) /2/ap]}? + {h + (v0/a)}* = 1/a%p?,4 


in which # is the parameter, the significance of which is that sin—(pvo) 
is the angle of emergence of the wave path at the surface. In fact, 


1 Published by permission of the Board of Directors, Humble Oil and Refining 
Company. 

2 Mathematician, Geophysics Department, Humble Oil and Refining Company, 
Houston, Texas. 

3 Cf., On Seismic Computations, with Applications, I, GEopHysics, pp. 9-22, 
Vol. I, No. 1, (1936). 

4 Loc. cit., equ. (12). 
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if, at any point along a wave-path we designate by 6 the angle between 
the positive 4-axis and the tangent direction of the path in the direc- 
tion of propagation, then 6=sin—(pv), in which v is the velocity at 
the point. The coordinate system used is that shown in Fig. 1, where 
the “‘well” is represented by the line «=d. 
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Perhaps the best method of obtaining the desired results in this 
case is to introduce a polar coordinate system for each wave path in 
the following manner:® 
(2) cos 6/ap = [(1 — po?)'/?/ap] — x, 

2 
sin 0/ap = h + (v/a). 
The travel time to the point P:(d,k) from 0 is given by 


0=0, 
(3) t= J ds/v, 


=; 


5 Here, as also in the next case, the analytic methods of the previous paper may be 
used directly. However, we felt that it is of some advantage to introduce the param- 
eter 6, both for purposes of making the integrations simple and to indicate other 
possible points of view. 
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in which 

(4) ds = (1/ap)d0 

and 

(5) v= %+ ah = sin 0/p. 


Accordingly, equation (3) becomes 
0» 
6) t= (1/2) q csc 6d6 
= (1/a) {log [sin 62/(1 + cos 62) | — log [sin 6:/(1 + cos 6) ]}. 
Obviously, 
(7) sin 0; = pv, sin 02 = p(vo + ah), 
in which / is the ordinate of the point P. However, we also have 
Cos 0, = (1 — pu?)1/2, 
cos 6. = F [1 — p2(vo + ah)?]1/2, 


where the negative sign is to be used for 62>7/2 and the positive sign 
for 0.<2/2;° ie., when h<[d?+(v,/a)?|"/?— (v/a), the negative sign 
is to be used, and, when + is greater than this quantity, the positive 
sign holds. The case when & is equal to this quantity corresponds to 
cos #.=0. Thus, equation (6) can now be written in the form 


t= (1/a) { cosh—"(1/ poo) + cosh! [1/p(v9 + ah) |} 
(9) = (1/a) cosh! {1 + (1 — p00?)/[1 — p2(v9 + ah)?]/?} 
/{ p?v0(v0 a ah)} , 


From equation (1), we find that the value of p corresponding to the 
path through P:(d, h) is: 


(10) p = 2d/{a?(d? + h2)? + 4v0(vo0 + ah)(d? + h®)} 1”, 
and, consequently, 
[xz — p2(v9 + ah)?]!/2 = F [ah? + 2hvo — ad?| 
/{a(d? + h?)? + 4v0(vo + ah)(d? + h?) } 1/2, 
Thus, finally, equation (9) becomes: 
(12) t = (1/a) cosh! {1 + [a2(d? + h?)/200(v9 + ah)]}. 


6 In any case, 0<@.<7. 


(8) 


(11) 
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Here we introduce the “dimensionless” quantities 
(13) r= at, o = ad/v, n = ah/v, 
and write (12) in the form 


(14) cosh r = 1 + (0? + n?)/2(1 + 7). 


The appearance of this family of curves is shown in Fig. 2, in which 
we have also inserted (in dashes) the locus of ‘‘minimum-time”’ 
points, whose equation is coshr=1-+7. 

o ; 4 4 4: s#s 2 4» 4 





FIG. 2 
THE EXPONENTIAL CASE: v= wexp(ah).? 
In this case, the family of wave paths has as its equation 
(15) sin! [ pv» exp (ahk)] — ax = sin~! (pv) .8 


7 The symbol exp(x) is used to indicate e. 
8 Loc. cit., equ. (26). 
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It is to be noticed that, in the equation, the angle 6=sin—[puexp(ah) | 
is measured in the same way as explained in the linear case. The value 
of the parameter p for the curve of the family passing through 
P:(d,h), (cf. Fig. 3), is 








(16) = sin ad/vo[1 + exp (2ah) — 2 exp (ah) cos ad]!/?. 
Accordingly, we use equation (3), in which, since 
(17) sin 0 = pvp exp (ah), 
(18) ds = (1/a) csc 6d0, 
and 
(19) v = (1/p) sin 8, 
O 
Gg, 











Oo 
\h v=veu 
FIG. 3 


we have 
6 


t = (p/a) “esc? 0d8 
(20) 6, 


= (p/a)(cot 6; — cos 62). 
Here, of course, 
cos 0, = (1 — pv?)!/?/ po 
(21) 
cot 02 = F [1 — pv? exp (2ah) |'/2/pu9 exp (ah). 
From (16), we have 
[x — pv? exp (2ah) ]!/2 = = [exp (ah) cos ad — 1] 


(22) /|x + exp (2ah) — 2 exp (ah) cos ad |1/2, 











304 M. M. SLOTNICK 


Consequently, (20) becomes 

(23) ¢ = (1/av) -_ (— ah)[t + exp (2ah) — 2 exp (ah) cos ad |1/?} , 
The ‘‘dimensionless” quantities that may be introduced here are: 

(24) T = avo, o = ad, n = exp (ah); 

and, as a result, equation (23) becomes: 


(25) t = (1/n)(1 + 9? — 20 cos o)'!?, 


These curves are shown in Fig. 4. It is to be noted that all the curves 
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have the line 7=1 as a common asymptote; that each of the curves, 
for which o<7z/2, has a minimum value for 7, and that the locus of 
these minima is the curve r=(1/n) (n?—1)"/?; that the curves for 
which ¢>z/2 have no minima but are rather such that 7 decreases 
monotonically with increasing 7. 

As a final point, we wish to remark that, considering ¢ as a function 
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of h, in well-shooting from a fixed point on the surface, the value of 
the derivative of ¢ with respect to / is given by the relation 
(26) dt/dh = cos 62/v. 


In this result, 62 is the angle between the positive h-axis and the direc- 
tion of propogation of the ray at the point in the well at depth 4, 
and v is the seismic velocity there. 








THE GROWTH OF COMPANY OWNED OPERATIONS IN 
GULF COAST GEOPHYSICAL PROSPECTING 
SINCE JULY 1930 


E. E. ROSAIRE! anp K. RANSONE! 


ABSTRACT 


Graphical data are presented for the refraction seismograph, torsion balance, and 
reflection seismograph in the Gulf Coast over the period July 1930 to February 1936. 
For comparison, the variations in the price of 30° Baumé crude oil at Houston are pre- 
sented for the same period. 

The results indicate that: (1) The refraction seismograph was discarded before the 
1931 minimum in the price of oil. (2) During this 1931 minimum, exploration of any 
kind was practically non-existent. (3) The 1933 minimum, on the contrary, seems to 
have effected contract operations more than company owned operations. (4) Company 
owned torsion balance operations seem to have reached a saturation point while con- 
tract torsion balance operations are rapidly decreasing from their 1934 peak. (5) Com- 
pany owned reflection operations have shown a steady increase of ten crews per year 
while contract operations have shown a very irregular increase, averaging six crews per 
year, since 1932. (6) New contract reflection companies have appeared since 1930 at 
the rate of three per year. (7) The peak of contract reflection operations will probably 
be reached in 1937 or 1938, after which these operations will disappear in the same 
fashion as did the contract torsion balance operations. 


INTRODUCTION 


Most of us are aware of the fact that company owned geophysical 
explorations have been increasing, and it is generally assumed that 
at some future time the geophysical contracting companies will dis- 
appear from the field. The available data have been analyzed here 
with a view to investigating this trend for one particular petroleum 
producing province, where the activity is greater than in any other 
part of the U.S.A., and where petroleum production is dominated by 
the operations of large oil companies like the Humble, Shell, Gulf, 
Texas, Magnolia and Stanolind. 

The data have been compiled from available scout reports. Since 
the making of these reports is a minor specialty of the local scout 
organization, the data are subject to frequent lapses, occasional glar- 
ing errors, and so are often erratic. As a result the charts show irregu- 
larities, generally due to the failure to report crews at all, or to re- 
porting a given crew operating at several places in the same week. 
Careful checking by the junior writer has eliminated many irregu- 
larities which appeared in the preliminary analysis, but the results are 
still far from being precise. This is to be expected, however, since the 
data are gathered from several sources, few of which are exact. 


1 Independent Exploration Company, Houston, Texas. 
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Further, deliberate attempts are frequently made to lead the scouts 
astray, and since the reporting of geophysical activities is only one 
field function of the company scouts in this area, they are to be con- 
gratulated that the data are no worse than they appear to be. In 
checking some hundred geophysical crews each week, mistakes are 
bound to occur. As a result we wish to thank Mr. Elton Rhine of the 
Empire Oil and Gas Company at Houston, Texas, for the information 
and reports that he has furnished us, and to congratulate him and 
his fellow scouts on doing a pretty good job in the face of material 
difficulties. 
THE REFRACTION SEISMOGRAPH 

Reports were available back to July 1930, so that we were able to 
get the fade out of the refraction technique. After plotting the price 
of oil over the period in question, we were justified in agreeing with 
the impression previously held by the senior author, viz., that the 
refraction technique was dropped because the law of diminishing re- 
turns ran its full course, and not that the 1931 decrease in the price 
of oil delivered the death blow. Just as Hobbs was not drilled on 
magnetic data because of the adverse results from the drilling of the 
three preceding magnetometer prospects, one outstanding refraction 
prospect on the Gulf Coast, later proved to be a large oil field by a 
wildcatter, was probably never drilled because of the three dry holes 
drilled just previously on three other refraction prospects. East Texas, 
however, effectively administered the coup de grace, and the dying 
gladiator never recovered. While the method was a new broom, how- 
ever, it swept clean and may have passed up only one shallow or 
moderately deep salt dome. Once the method had shot its wad, how 
ever, it was through, and East Texas helped make that fact apparent. 


THE TORSION BALANCE 


We regret that data back to 1924 or 1925 are not available on 
torsion balance activity, for the curve back of 1930 would be most 
interesting. It would show a low corresponding to the point where 
the gravity maxima finally proved disappointing, and when the im- 
portance of gravity minima had not yet been appreciated. 

From the data at hand, however, it is of interest to notice how 
the decrease in total torsion balance activities from late 1930 to the 
middle of 1931 follows the drop in the price of oil during the same 
period, but after that date an effect due to the price of oil does not 
seem apparent, for the 1933 drop in the price of oil is not reflected on 
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the torsion balance activity chart. However, at the first of 1933 a 
change from zero increase to a definite increase in company owned 
torsion balance activity seems to coincide with the initial production 
at Livingston, a gravity prospect. Paradoxically, Livingston was al- 
most the last of the unassisted gravity discoveries, being followed only 
by Garwood and Cleveland. Production at Cleveland in late 1933 
seems to coincide with an increase in contract torsion balance ac- 
tivity, carrying up to the peak which was reached and passed by late 
1934. 

About the time of the discovery of production at Hastings, the 
recognition became prevalent that torsion balance prospecting pointed 
the way for the seismograph, not the drill. As a result, the conclusion 
was drawn that isolated torsion balance surveys were of little or no 
practical value, and by common consent the bulk of torsion balance 
exploration has been left to the larger operators who can see their 
ways clear to mapping county by county, thereby turning up sub- 
marginal gravity prospects for later evaluation by the seismograph. 
As a result, contract torsion balance activity has reached the vanish- 
ing point, for although high grade equipment has been on the market 
for some time, many of the men who gained material torsion balance 
experience as contractors have been drafted by the large operators 
of their own crews. 

Any revival of torsion balance contracting is practically out of the 
question, for those companies who do not have a material background 
of experience and data are hardly justified, in fact, there is not time 
enough remaining, to duplicate the extensive net works already com- 
piled by the Humble, Shell, Gulf, Magnolia, Texas and Stanolind. 
However, the availability of a fieldworthy gravity meter might tempt 
some of the other operators into an extensive program, for the rate 
of coverage would be materially increased thereby, at a reduction in 
cost which would probably be considerable. 


THE REFLECTION SEISMOGRAPH 


Although the reflection seismograph had already been introduced 
into this province in late 1930, the price of oil was falling, and ex- 
ploration was decreasing in popularity. The initial success of the 
method followed its use in the detail of the Iowa prospect, then recog- 
nized as a torsion balance minimum, and as something better than a 
marginal refraction prospect. The first well drilled was located on the 
reflection apex, and fortunately produced, for there seems to be no 
well developed central graben. 
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Other successes were largely in the condemning of prospects orig- 
inated by other methods, and several dry holes confirmed early pre- 
dictions of the reflection seismograph. Though there was little or no 
demand for such a method, still an excellent foundation was laid for 
its future use, for since the recognition of the possibility of mapping 
by dips had been made at Darrow, the idea was slowly but generally 
becoming prevalent that this was a workable method. 

When the importance of the Conroe field was realized, there was 
large immigration of Oklahoma oil operators into the coast. Many of 
them were familiar with the recent and spectacular successes there of 
the reflection method as introduced by the Geophysical Research 
Corporation. These newcomers did not show the sales resistance that 
many of the resident companies did, and thereby the introduction of 
the method was facilitated. Note, however, that some company 
owned reflection operations had begun in 1930, and credit must be 
given to the Humble for having appreciated the possibilities of the 
method at such an early date. 

After Conroe, reflection methods were generally adopted, and the 
total operations increased at an average rate of fifteen crews per year 
until at present there are some seventy-five crews, all told, operating 
in this province. 

It is interesting to note that the growth of company operations 
was not affected by the 1933 price cut, while the contract activities 
show nearly a fifty percent decrease at the same time. Note also that 
the company owned operations show a remarkably steady increase 
each year, while the contract operations are nowhere near as regular. 
Further, the average (total) annual increase is made up of a company 
owned increase of not quite ten per year, leaving a contract increase 
of less than seven per year, but this is by no means a steady increase. 

The material increase of total operations in late 1935 occurred 
largely in Chambers, Jefferson and Orange Counties (see regional 
activity report for Southeast Texas and South Louisiana, Geophysics, 
Vol. I, No. 1, Jan., 1936) and seems to have been due to the discovery 
of oil at North Anahuac, and the recognition of the importance of 
faulting. There is some suggestion that the usual fall off, after this 
peak, was eliminated by the crude oil price increase at the start of 
1936. 

The usual exploration history of any area shows a distinct over 
run of prospecting, a final phase of hysterical exploration running 
past the peak of discovery. California is an excellent example of this 
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determination to exhaust every possibility of finding structure. How- 
ever, there was not time enough there to build up appreciable com- 
pany owned operations before this over run took place, so that there 
the contractors enjoyed a period of great activity. 

However, as the chart shows for the Gulf Coast, company owned 
operations in reflection exploration were initiated shortly after the 
introduction of the method by the contractors, and have grown at a 
rate greater than have the contract operations. Of the eleven largest 
users of the method on the Gulf Coast during 1935, the two largest 
were company owned. As a matter of fact, of the nine largest users, 
only one entirely, and only two partially, relied upon contract opera- 
tions. There has been a distinct trend towards company operation, 
the last major defection from the contract using ranks occurring 
when the Stanolind bought the Western Geophysical Corporation a 
year or two back. 

In view of this, and with the torsion balance statistics in mind, it 
is quite probable that the reflection seismograph over run in the Gulf 
Coast will be made largely by company owned crews. At the present 
rate it is reasonable to predict that 1937 or 1938 will be the peak year 
for contract operations in this province. 

In conclusion, it might be well to examine the validity of this 
trend towards company owned operations. In general, it is based upon 
the following ideas. 

1. By having the work fully organized within the limits of one 
organization, smoother control will be effected, and therefore, less 
opportunities for leakage of important information. 

2. The contract companies, more or less openly, are interested in 
finding oil fields of their own, thereby competing with their clients. 

In regard to the first reason, we point out that a large company 
organization is no bar to leaks, and are of the opinion that, on the 
whole, the leaks within company organizations have added up to ma- 
terially more embarrassment than have those leaks within contract 
organizations. In regard to the second, we point out that, as serious 
competitors in this game of finding oil fields, the total interest plays 
of all the contracting companies do not result in anywhere near the 
competition furnished by any one of the larger users, Shell, Humble, 
Gulf, etc. Under the circumstances, the second reason, although one 
which has been advanced in all seriousness, is much the case of chok- 
ing on a gnat and swallowing a camel. 

Of course, adding a group of one, two or three reflection crews to 
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an oil company is, apparently, a relatively simple job these days. 
Equipment capable of recording reflections has been on the market 
for some time, and the contracting companies have many men, ca- 
pable of operating this equipment, who could be made available at the 
necessary salaries. In fact, this idea has been carried to the point 
where equipment is available at very nominal rentals without any 
restrictions on the interpretation of the data secured by observers 
who record under specific instructions, and do nothing more. 

However, there is something to be said for many contracting com- 
panies. Some of them go to extremes in order to turn out work of high 
caliber. The worst that can be said of these companies is that they may 
take an over optimistic attitude in their desire to avoid passing up oil 
fields. At that, perhaps, the dry hole record of some of these con- 
tracting companies is no worse than that of company owned organi- 
zations, and sometimes their record of discoveries is more impressive. 
These contractors, in general, do a volume of business which is large 
enough to justify the maintenance of appreciable and frequently 
capable research staffs, a function which, as Dr. Weatherby pointed 
out in his paper on “The Organization of an Effective Exploration 
Department’”, they should carry out. 

Dr. Weatherby also pointed out how a contract organization ties 
into a properly organized exploration department. It is frequently 
true that some contract companies are better organized, and often 
better staffed than some of the competing company owned organiza- 
tions. Therefore, our suggestion is that before the executives of any 
company act on a decision to tack on a few reflection crews and so 
drop the contractors, they should examine carefully the existing 
organization of their own exploration department, and then consider 
again, that in view of the continuing progress in the art, there is 
probably a minimum number of reflection crews which can be oper- 
ated in order to justify the organization and maintenance of the re- 
quired research and auxiliary personnel. 


2 Geophysics, Vol. 1, No. 2, June 1936, pp. 179-188. 











EFFECT OF DEVELOPMENT TIME AND DEVELOPER 
TEMPERATURE ON THE PRODUCTION OF PHO- 
TOGRAPHIC SEISMOGRAPH RECORDS 


FREDERICK A. TOMPKINS! 


ABSTRACT 


The effect of variation of development time and developer temperature on the 
production of uniform records is discussed, and a method of determining the relation 
is shown. Graphs plotted from sensitometric data indicate the need of some means of 
temperature control. 

In a previous paper by the author,? some of the fundamental 
factors influencing the photographic processing of seismograph records 
were briefly outlined. Following that presentation, the desire was 
expressed for a more detailed and graphical presentation of some of 
these relations. In this paper the most important relation of develop- 
ment time and developer temperature to effective contrast will be 
discussed. 

This relation can be determined most accurately by sensitometric 
methods. If a strip of photographic paper be exposed to a constant 
light source for a series of increasing exposure times, a graded series 
of silver deposits is secured’ on development. The density of these 
deposits can readily be measured, and the values can be expressed 
mathematically by the relation 


D=log 1/R 


where D=density or blackness of the deposit and R=per cent of the 
incident light reflected from the silver deposit. On plotting these 
density values against the logarithm of the corresponding exposure 
values, the typical characteristic curve of a photographic material is 
secured (Fig. 1). 

Now by developing several identically exposed strips for increas- 
ing times at a given developer temperature, and then by plotting the 
measured densities, a family of curves, Fig. 2, is secured. Such a 
family of curves shows the rate. of density increase as a function of 
development time. 

It will be noticed that the maximum density of this sensitive ma- 


1 Paper Manufacturing Division—Eastman Kodak Co. 
2 Fundamental Photographic Processing Operations Influencing Production of 
Seismograph Records. Geophysics, Vol. 1, No. 1, Page 107, Jan. 1936. 
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Fic. 1. Typical curve showing the relation between density and log exposure. 
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terial is not attained until the development time approaches 60 
seconds when development is carried out at 50°F. At 70°F. (Fig. 3) 
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Fic. 3. Compare with Fig. 2 to note effect of raising developer temperature. 
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Fic. 4. Compare with Figs. 2 and 3. Note appearance of fog at 60 sec. development 
as indicated by the lower part of the curve lifting from the zero density line. 


the maximum density is approached at 30 seconds, and at go°F. 
(Fig. 4) this value is below 20 seconds. But, if the development time 
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is prolonged to 60 seconds, an over all greyish deposit or fog is pro- 
duced. With the developer at 110°F., development is completed under 
ten seconds, for at this time the maximum density has been obtained, 
and the fog value has become pronounced. (Fig 5). Note how rapidly 
the fog density increases with increased development at this tem- 
perature. 

Hence, it will immediately be seen that such an increase of fog 
produces a degraded image quality. Accordingly, if we subtract the 
minimum density from the maximum density (D max.—D min.) for 
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Fic. 5. Compare with Figs. 2, 3 and 4. Observe the rapid increase of fog 
as development time is increased at this high developer temperature. 


each curve of these several families, we secure a numerical value of 
the effective contrast of the image in each case. The term, effective 
contrast, is used to denote the optimum condition of blackest black 
with the whitest white obtainable under each development condition. 
In short, it is a numerical valuation of the best image that the 
operator can produce under a given set of conditions. 

Again, if we plot these values of effective contrast as a function 
of development time for each developer temperature, a series of 
curves (Fig. 6) is obtained which indicates the rate of change of 
effective contrast and which begins to show the relation of time and 
temperature to the final result. 
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Carrying the reasoning one step farther, if the development time 
is determined from the curves of Fig. 6, for any constant value of 





CHANGE OF CONTRAST WITH TIME & TEMPERATURE 
14} INSURANCE BROMIDE R D-72 DEVELOPER I-I 





EFFECTIVE CONTRAST 
é @ 
' 














1 lL 1 l l i i] | 
10 20 30 4Q@ 50 60 70~=—s 80 90 
DEVELOPMENT TIME - SECONDS 
Fic. 6 


effective contrast, and if these development times are plotted against 
the developer temperature, the final time-temperature relation be- 
comes known. Fig. 7 shows the relation for Insurance Bromide Grade 
R, and Fig. 8 depicts the relation for Recording 809 when developed 
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in D-72 diluted 1-1. Both papers are commonly used in seismograph 
recording work. 

Both curves indicate the great need for adjustment of develop- 
ment time in order to produce constant results over a wide range of 
temperature—a range not uncommonly found in actual field work. 
More important, however, they indicate the need for some means of 
temperature control, for it becomes practically impossible to control 
development at the higher temperatures, while the prolonged time 
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necessary to produce satisfactory images at the lower temperatures 
may retard field operations unnecessarily. 

The relations shown are for fresh developer. Used developer shifts 
the curve vertically along the time axis, at the same time increasing 
its steepness in proportion to the degree of exhaustion. 


SUMMARY 


The method of determining the development time-developer tem- 
perature relationship for two recording papers has been discussed, 
and graphs indicating such relationships are shown. These graphs in- 
dicate the need for temperature control. It is suggested that these 
curves can be used as the basis for determining actual operating pro- 
cedure. 














THE AMPLITUDE AND CHARACTER OF 
REFRACTION WAVES 


A. WOLF* 


ABSTRACT 


In seismic exploration it is assumed that elastic waves can be refracted into and 
out of a high speed medium at the angle of total reflection and meanwhile travel along 
the interface between the two media. That such a path is in agreement with the general 
theory of elastic waves has been proved by Muskat. The present paper is an extension 
of the work of Muskat, and shows in addition that the general theory can explain the 
low amplitude and low frequencies of these refracted waves, and is otherwise in agree- 
ment with observation. 

Equations are developed for the effect at any point of the low speed medium of an 
impulse or of a periodic disturbance initiated at any other point of the same medium. 
Details are worked out for the case where the recorder is so situated that the refracted 
wave arrives first. The direction of the arriving wave is found to be that of the angle 
of total reflection, and its amplitude to diminish with the square of the distance from 
the shot point—in agreement with Muskat. The amplitude of the refracted wave is also 
found to be only a small fraction of that of the direct wave, or of that of the reflected 
wave which arrives still later. The motion produced by the refracted wave is found to 
be a unidirectional pulse which lasts until the arrival of the direct wave. 

Finally, the appearance of a record of motion is given for certain assumed numerical 
data. For this case it is found that the average amplitude of the refracted wave is 
8 per cent of that of the direct wave. 


It is generally assumed in seismic exploration that elastic waves 
can berefracted bothinto and outof a high speed medium at the angle 


LOW SPEED MEDIUM 
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of total reflection, to be recorded on the surface of the ground. A 
path, such as shown in Fig. 1, has become associated with this work. 
The object of the present paper is to show that the general theory of 
elastic waves accounts for the existence of that type of propagation 
and is in general agreement with observation. It will be possible to 
explain both the low amplitude and the relatively low frequencies of 
these refracted waves. 

Perhaps the first paper to call attention to the application of the 


* Geophysical Research Corporation. 
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exact equations to this subject was published by Muskat.' He used a 
method well known from the work of Sommerfeld on the propagation 
of electromagnetic waves. Muskat’s work proves that the type of 
path assumed in seismology for refraction waves really exists, but his 
conclusions about the order of magnitude of the amplitude of such 
waves seem to be incorrect, and his discussion gives very little idea 
of their character. 

The present discussion follows the methods developed in a paper of 
H. Lamb: “On the Propagation of Tremors over the Surface of an 
Elastic Solid.’” 

Assume that there are two elastic solids A and B extending to in- 
finity in all directions from a plane of separation SS between them 
(cf. Fig. 2). The reciprocals of the velocities of propagation of 


4° 





OR 
j 





L_ ny S-y 





Fic. 2 


longitudinal waves in A and B will be denoted by a and 6 respectively. 
It is assumed that a is greater than b, which means that the velocity 
in A is smaller. Writing 

a/b=n 
it follows m>t1. 

Let the densities of A and B be respectively y4 and yz. A system of 
cylindrical coordinates (7, z) will be used whose origin, O, is located on 
SS. The positive direction of the z axis is in A. At a distance h above 
O there is a small cavity O’ in which a shot, which is supposed to be 
radially symmetrical, is set off at the time ¢=O. At a point Po, a 
distance R, away from O’, there is observed the elastic displacement 
u(t). According to the hypothesis made regarding the nature of the 
shot, this displacement is directed along the line O’P,. In what fol- 
lows the function u,(¢) will be referred to as the impulse. The prob- 


1 Muskat: Physics 4, 1933, Pp. 14. 
2 Lamb: Phil. Trans. Roy. Soc. 203, 1904, p. 31. 
3 The effects produced by the surface of the ground are not considered in this paper. 
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lem is to determine the horizontal and vertical elastic displacements 
v, w, at a point P (coordinates 7, z) in A, which correspond to the 
impulse u,(¢). It is somewhat simpler to study the refraction as a 
function of u,(¢) than it would be to derive the corresponding equa- 
tions if the pressure in the cavity O’ were given as a function of time. 

The procedure can be outlined as follows: First the expression for 
the elastic potential of longitudinal waves ¢ is determined on the 
assumption that the impulse is a periodic function of time 

u(t) = Uper® (tak o) | 

By the use of the Fourier integral an equation for ¢ is then obtained 
which corresponds to an impulse of infinitesimal duration (cf. Fig. 3). 
By superposing such solutions it is easy to determine ¢, v, w for any 
impulse u,(t). 











dT 
+6 
T — 
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FIG. 3 


In the case of a periodic impulse, it is easy to show that 


uoRo (* Adr - 
o¢ = — ———— fo f dx. 
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(1) nr J_,VN-—1 
exp jw[t — V1 — A®a(z + h) — Aarcosh x] 
where f (A) is the coefficient of reflection of plane longitudinal waves 
at an angle ¢ given by the equation 
sin g=A 


I t 
I +f dt. 
aRod o 


Except for minor changes, Equation (1) is the same as Muskat’s 
equation for the potential. The Bessel function which he uses is re- 
placed by an integration with respect to « in Equation (1). The 
method of deriving such equations is well known from Lamb’s work 
on Rayleigh waves and will, therefore, be omitted. The function 


and a is the operator 
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f(A) depends on the Poisson’s ratio of the solids and is very compli- 
cated. If, however, the transformation of longitudinal waves into 
transverse is neglected, the expression 


wis — i ~ yale — iP 
ye? — 1 + ya? — 1/2? 





(2) fA) = 


is obtained. In what follows f(A) will be assumed to have this simple 
form and furthermore ys =vyz. It is possible to calculate numerically 
the coefficients of reflection is special cases, and show that no great 
differences exist between the exact values of the reflection coefficient 
and form given by Equation (2), at least in the range of values of \ 
which is relevant. The following equations could, to a great extent, 
be written out with f(A) undetermined, and only the final calculation 
requires its knowledge. This makes it possible to state that no qualita- 
tive differences exist between the conclusions reached in this paper 
and anything which could be obtained by using more exact expres- 
sions. 

The integrand in Equation (1) can be regarded as a function of 
the complex variable \, provided, only, care is taken to use a suitable 
path of integration in the \-plane. Such a path is shown in Fig. 4a. 
The square roots must be taken positive for \=+ . Since the inte- 
grand vanishes exponentially at infinity in the lower half-plane, it is 
possible to modify the path of integration of Fig. 4a into that shown 
in Fig. 4b, which is chosen so as to make all subsequent integrations 
convergent. It will be noticed that the point 
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— Vr? + (zg + h)? 








is shown between 1/z and 1. This means that the recording station P, 
is located in such a way that the refracted wave reaches it before the 
direct, which is obviously the only interesting case. If transverse 
waves were taken into consideration, the contour of integration would 
have to be modified only slightly. The most important change would 
be the addition of small circles around the poles of f(A) which cor- 
respond to waves similar to Rayleigh Waves propagated along the 
interface of A and B. 

The refracted waves are given by those parts of the contour of 
integration which are between 1/z and a. The potential is given by 
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oo, = — m a. 
(3) = 


exp jw[t — V1 — Fate + h) — ar cosh x]. 





Multiplying the right hand side by dw and integrating from O tom, 
the potential corresponding to an impulse of short duration dT and 
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amplitide u, (cf. Fig. 3) is obtained. Since the quantities between 
the square brackets are real, the result of an integration of the ex- 
ponential with respect to w is 


5[t — V/1 — A2a(z + h) — dar cosh «| 

where the function 6 is defined by the equations 
dT \? 
i(y) =1 for y< (—) 

2 


4 The half circles around 1/n and 1 in Fig. 4a should be on the upper side of the 
real axis. 
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dT \? 
5(y) =o for y?> (—). 


2 


After this has been done, it is easy to carry out the integration with 
respect to x and obtain a relatively simple equation for the elastic 
potential. Differentiating it with respect to 7 and z gives the elastic 
displacements 
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and v=w=o for t<é,. 
\. is a parameter such that 


1/n<do <1 and t= V1 — Agrta(z + A) + Acar 


t, is obviously the time at which the refracted wave reaches the record- 
ing point P. It is, of course, in agreement with calculations made on 
the basis of the minimum time principle. 

The evaluation of the integrals is much simplified by the assump- 
tion 

z+th=o 

which means that both the shot point and the recording apparatus 
are located close to the surface of separation of A and B. The results 
must, apart from details, be true when (a+h)/r< <1. 

The following expressions are obtained in this case 
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I I 
for ¢>—ar and v=w=o for t<—ar 
nN nN 
where 
t?/a®r? — 1/n? 


k= 
1 — 1/n? 





and F(k) and E(k) are respectively the complete elliptic integrals of 
the first and second kind. It follows that, at the time of arrival of 
the wave, the direction of motion makes an angle with the normal to 
the surface of separation SS between the two solids, which is equal 
to the angle of total reflection. 

It is seen also that the amplitude of motion decreases as the 
square of the distance from the shot point; this is in agreement with 
the results obtained by Muskat. 


I 


aRop 





It remains to explain the nature of the operator o=1+ 


and to generalize the results for the case of pulses of finite duration. 
First, as regards the operator, it is possible to show that, if P, is not 
too close to the shot point and the time is not much greater than aR,, 
the operator can be neglected. For any other case, the operator can 
be evaluated according to rules of the operational calculus. 

It is very easy to generalize Equation (5) for the case of impulses 
of finite duration. They are first rewritten in the form 


Ro 
v= — f(t)udT 
ar? 


Ro 
w = — g(t)uodT.® 
ar? 


If the actual displacement at R, is u,(¢) instead of u,6(¢—aR,), 


Ro (* 
v= =| uo(x) f(t — x)dx 


Ro 6 
w =f uo(x)g(t — x)dx. 
ar? 0 


Let it next be supposed that the impulse is confined to a short, 
though not infinitesimal, time 7,. Then, approximately 


5 f, g are of the order of magnitude of 1. 
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Ro To 
v= = iy f Uo(x)dx 


Ro To 
w= — co f uo(x«)dx. 
ar? 0 


The character of the motion produced by the refracted wave is seen 
to be independent of the shape of the impulse and consist of a long 
pulse occupying the whole interval of time between the first arrival 
of the refracted wave and the arrival of the direct wave. Aside from 
a numerical factor of the order unity, the ratio of the amplitudes of 
elastic displacement of the refracted and direct waves is given by 
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FIG. 5 


In conclusion, the appearance of a record of motion will be dis- 
cussed for a case in which 7,=0.02 sec., m=1.5, 1/a= 10,000 ft/sec., 
r=5o000 ft. A sketch of the record of horizontal motion is shown in 
Fig. 5. The refracted wave arrives at the recording station at the 
time 0.33 sec. after the shot and consists of a unidirectional pulse 
which lasts until the time 0.5 sec. when the direct wave arrives. Its 
average amplitude is 0.08 the amplitude of the direct wave. A short 
time after the arrival of the direct wave a reflection from the interface 
is observed. Since the angle of incidence is assumed to be very nearly 
grazing, the direction of motion of the reflection is negative. 














THE SEISMIC ELECTRIC EFFECT! 
R. R. THOMPSON? 


ABSTRACT 

The seismic electric effect is the name which has been given to the variation of earth 
resistivity with elastic deformation. An account is given of the use of this effect in 
recording seismic waves. The recording circuit employs a direct current sent through 
the ground by means of a storage battery and spaced electrodes, and the changes in 
current due to seismic disturbances are recorded. An analysis of the noise appearing 
in the recording circuit is given and methods of reducing the noise level are described. 
The results of measurements on the magnitude of the effect are given, and the sensi- 
tivity obtainable with the seismic electric method is compared with that given by an 
electromagnetic pickup. 

The possibility that there might exist a variation. of earth resistiv- 
ity with elastic deformation was first suggested by Louis Statham and 
L. W. Blau of the Humble Company. They pointed out that, if the 
effect existed, it might be used to provide an entirely new method of 
picking up the seismic energy in reflection and refraction shooting. By 
a suitable arrangement of electrodes it would be possible to obtain a 
current distribution in the earth involving almost any desired volume 
of earth. Since the resistivity at any point in this volume would vary, 
as a function of time, with the elastic strain at that point, the varia- 
tion in effective impedance between the two groups of electrodes 
would be the integral, throughout the volume, of the seismic dis- 
turbance as a function of time. Due to the very high apparent velocity 
(along the surface) of the reflected waves as compared to waves of 
less penetration, the method would offer a possibility of discriminating 
against the shallow waves and in favor of the reflections. 

This effect, the variation of the earth resistivity with elastic de- 
formation, which has been termed the seismic-electric effect, has been 
investigated as to causes, magnitude and possibilities of use. 

A schematic diagram of the circuit used in recording the seismic 
electric effect is shown in Fig. 1A. For simplicity, only two electrodes 
A and B are shown, although either of these may be a multiple elec- 
trode earth contact. A direct current is passed through the earth be- 
tween these electrodes by means of a storage battery in series with 
the primary of a transformer. The secondary of the transformer is 
connected through an amplifier to an oscillograph and the signals are 


1 Published by permission of the Board of Directors, Humble Oil & Refining 
Company. 
2 Geophysics Department, Humble Oil & Refining Company. 
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recorded as in ordinary reflection shooting. Fig. 1B represents the 
equivalent input circuit in which R represents the effective ground 
resistance, Z is the impedance looking into the transformer primary, 
and J is the direct current flowing in the circuit. The resistance R is 
assumed to vary with time as 


R = Ry + Ar(t) 
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where Ar(¢) is very small compared to Ro. This circuit is approxi- 
mately equivalent to that of Fig. 1C where the effect of the variation 
in resistance is represented by a series generator of voltage Ar(t)I 
and internal impedance Ro. Thus the signal voltage effective in the 
input circuit is given by the product of the instantaneous value of the 
resistance variation and the magnitude of the direct current flowing 
through the ground. 

Noise voltages which originate in the input circuit may be divided 
into two classes: Noises due to resistance changes other than the 
seismic electric effect, such as variations in contact resistance; and 
noises due to currents already flowing in the earth which produce a 
fluctuating potential difference between the electrodes. The first type 














THE SEISMIC ELECTRIC EFFECT 329 


of noise increases linearly with the ground current so long as the 
current is not raised to sucha high value that the unwanted resistance 
fluctuations cease to be independent of the current. The second type 
of noise voltage is independent of the ground current and is a function 
only of the location, separation, etc., of the electrodes. In addition, 
it has been found that for large values of current a third noise voltage 
component must be included. This component is negligible for small 
values of current but increases rapidly as the current is raised above 
a certain value. Apparently this source of noise is associated with 
chemical action at the electrode. 
Thus an expression for the noise voltage may be written 


E, = IR, + E, + G(D) 


where 7 is the ground current, R, represents resistance fluctuations 
other than those due to the seismic electric effect, E, represents the 
noise voltage picked up due to natural or artificial earth currents, 
and G(J) is a term included to take care of the very rapid increase of 
noise with ground current when the current becomes too high. 

These noise components are illustrated in. Fig. 2 where it is seen 
that at very low currents EF, is the important component; as the 
current is increased JR, becomes the important component; and for 
very large values of current G(J) becomes the principal source of 
noise. Since the signal voltage varies linearly with the ground current, 
that is 

E, = Arl 
it is evident that the best signal to noise ratio is obtained in the region 
where the /R, term is the principal noise component and that in this 
region the ratio is 


The first experiments with electrodes showed that any kind of 
metal rods driven into or buried in the earth gave an intolerably high 
noise level due to fluctuating contact resistances. A much quieter 
earth contact was finally obtained by drilling holes in the ground with 
a hand auger, filling the holes with a salt solution to serve as an 
electrolyte, and placing the electrode in the solution. The electrodes 
used were cylinders of copper or zinc about 33 inches in diameter and 
3 to 6 inches in length. Electrodes of this type were used in all sub- 
sequent work. 
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Measurements were made, with two such electrodes, on noise as 
a function of electrode separation. The results are shown on Fig. 3. 
Line A represents the equivalent input noise voltage for various dis- 
tances between electrodes. It is evident that in this case the noise 
increases linearly with distance between electrodes. This indicates 
that the noise comes from ground currents which are distributed in 
such a way that the potential gradient is very nearly constant over 
the area involved. This would be expected to be the case if the earth 


FIG. 2. 


ILLUSTRATION OF NOISE COMPONENTS 
IN SEISMIC ELECTRIC RECORDER. 


EFFECTIVE INPUT NOISE VOLTAGE. 














GROUND CURRENT. 


is homogeneous, and the source of the disturbance is not a local 
source near the electrodes. If this is true, it should be possible to 
balance out the voltage picked up in this way by using the arrange- 
ment B, where three electrodes are set in a straight line, the center 
electrode being equidistant from the two outside electrodes. The two 
outside electrodes are connected together and are of one polarity 
while the center electrode is of the opposite polarity. With an arrange- 
ment of this type the noise dropped to the value indicated by ‘‘B”’ 
on Fig. 3 when the distance between outside electrodes was 50 feet. 
Five electrodes in line gave a noise voltage indicated by ‘“‘C” on 
Fig. 3. Any number of electrodes may be used in a balanced arrange- 
ment, it- being necessary only to space the electrodes so that the 
vector sum of all the lines correcting positive electrodes to negative 
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electrodes is zero. It thus appears that when only two electrodes are 
used the £, term is relatively large and increases linearly with elec- 
trode separation. With a balanced arrangement the E, term can be 
reduced considerably. 

The results of measurements on the variation of noise with ground 
current are shown on Fig. 4. In this case a three electrode balanced 
circuit was used so that the EZ, term should be small, and the 7R, 
term should be the important one. It will be noted ‘hat the noise 
increases linearly with the current and has a very small value at zero 
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current as would be expected, since the noise at zero current should 
be the value of £,. 
As has already been shown, the effective input signal voltage is 
given by 
E, = Arl. 


The value of Ar, and hence of the signal voltage, varies with time 
in exactly the same way as the ground motion. By inserting a low 
impedance generator of proper frequency and known voltage in the 
input circuit and comparing the amplitude recorded from this input 
to that obtained on the record of a dynamite shot, the magnitude of 
E, may be determined at various times on the record. The magnitude 
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of the ground current may be measured by a series ammeter and the 
ground resistance determined by the 7R drop across the electrodes. 
Thus Ar and R may be determined and the ratio Ar/R computed. 
This ratio is, of course, a function of time along the record, and for 
any particular set-up will vary with the size of the charge in the same 
way as the ground motion. Fig. 5 shows the values obtained for 
Ar/R for a number of locations in the Gulf Coast. The values of Ar 
are taken from the maximum amplitudes recorded in the neighbor- 
hood of .5, 1, 1.5, and 2 seconds. It will be seen that the resistance 
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variation dué to the energy arriving at about one second out on the 
record is of the order of 1 to 15 parts in 10 million and at two seconds 
it is about a tenth as much. 

The charges used vary from 2} to 10 pounds and are of about the 
same size as the average charge for ordinary reflection shooting in the 
same locality. The distance from the shot to the recorder varies from 
2000 to 2500 feet. The largest effect observed, that of No. 4, represents 
a usable amplitude on the record about equal to that obtained from 
a standard electromagnetic pickup. 

Photographs of two records with seismic electric traces are shown 
on Fig. 6. On record No. 1 a standard pickup was placed beside the 
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center electrode of the seismic electric set-up which was 1900 feet 
from the shot; the outside electrodes were fifteen feet on either side 
of the center electrode. Identical amplifiers of the resistance coupled 
type with no low frequency elimination filters were used for the two 
traces. On record No. 2 the fourth trace of a standard field set-up was 
replaced by a seismic electric trace. The spread for‘this record was 


SE x 107 





2100 to 2850 feet, number three pickup was at 2400 feet and the 
seismic electric electrodes were at 2375, 2400, and 2425 feet. The 
sensitivity on these records is sufficiently high for practical field use 
but they represent an effect of about the magnitude given by curve 4 
of Fig. 5, that is, they were taken in the most favorable areas found. 
In many areas a much lower sensitivity was found to exist. The 
problem of increasing the usable sensitivity is entirely a matter of 
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reducing the noise level and it is possible that means for accomplish- 
ing this may be found at a future date. 

Since the potential gradient, due to the current sent through the 
earth, is large near the electrodes and drops off rapidly for increasing 
distances from the electrodes, the principal contributions to the input 
signals come from the regions near the electrodes. Consequently the 
seismic electric recorder does not integrate the seismic disturbance 
uniformly throughout the volume of earth which carries the current, 
but gives most weight to the disturbance near the electrodes. 

As to the mechanism of the resistance change, no definite answer 
has been found. Experiments on electrodes set close together so that 
the intervening earth becomes saturated with the electrolyte show 
very little response, indicating that the effect is in the earth itself and 
not in the electrolyte or at the electrode surface. It has been sug- 
gested that the effect might be due to the small temperature fluctua- 
tions which occur with the compression in the earth and that, on 
account of the temperature coefficient of resistivity of the earth, a 
resistance fluctuation is produced. It seems unlikely, however, that 
if this were true, there would exist such large variations in the magni- 
tude of the resistance change as have been observed from one location 
to another. A more probable explanation seems to be that the effect 
is some sort of loose contact phenomenon between the particles of 
the earth. 











A SIMPLIFIED CIRCUIT OF THE SEISMIC ELECTRIC 
METHOD AND ITS STEADY-STATE SOLUTION! 


M. M. SLOTNICK? 


ABSTRACT 

The Seismic Electric Effect gives rise to the problem of finding the steady state 
response of a circuit consisting of an inductance and a response of a circuit consisting 
of an inductance and a resistance of the form R+A cos cot (R>4A) in series with a D.C. 
input. In this paper a solution is given, other than the one usually obtained by the 
method of successive approximations. 

The circuit shown in the accompanying figure was suggested by 
the seismic-electric method? as being the simplest possible one which 
might arise in a very ideal situation. Although the amplitude and 
phase of as many of the harmonics as desired can be computed to any 
degree of accuracy by the usual and well-known methods? of successive 
approximations, it may not be entirely without interest to indicate 
the solution in series form and the method by which this solution 


is obtained. 











Fe + Acoswt (R>A) 


In the figure, Z denotes a constant inductance in series with a 
d. c. input voltage EZ, and a resistance of the form R+ A cos wf, in 
which R and A are both constant with R>A. The differential equa- 
tion for the current J is 


1 Presented by permission of the Board of Directors of the Humble Oil and Refin- 
ing Company, Houston, Texas. 

* Mathematician, Geophysics Department, Humble Oil and Refining Company. 

3 Cf. R. R. Thompson, The Seismic-Electric Effect and its Use in Recording 
Seismic Waves, Geophysics, Vol. 1, No. 3, pp. 327-335. 

* Cf.; e.g., E. A. Guillerin, Communication Networks, vol. 1, (New York, 1931), 
Pp. 403, et seq. 
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(1) i Aes ee 
t 


and the solution, subject to the boundary condition, 


(2) I(o) =o 

is 
t 

(3) T= peratovines f erttbsinw td, 
0 


in which we have set 
(4) k= E/L, a= R/L, b = A/Lo. 


The problem, then, is to find J(¢) as >, inasmuch as we are 


interested in the steady-state part of (3). 
Let us assume that this steady-state part of J(¢) is 7*(#) and that 


its Fourier expansion is 


(5) I*(t) = >> (in cos nwt + jn sin not). 


n=0 
Since J*(#) must satisfy (1) identically in ¢, the Fourier coefficients 
i, and j, must satisfy the equations: 
aig + Zbwi, = hk, 
ai, ++ wt + bwi + A bwie = Og 
(6) = wt ot. ai -}- sbw/» =O, 
Qin + NOG n + 5 bw (int + in+1) == Os 
, 1 , ‘ mS) a Ge 
aJn — NW In, + 5bw(jn—1 + jnti) = Che 


If, then, we succeed in finding the series (5) for 7*(¢), we shall also 
succeed at the same time in finding the solution for this infinite set of 
linear equations in the infinite number of unknowns 2@y, jn. 

We begin with the well-known expansion®: 


(7) e~*sinet = Jy + 2). Jon cos 2mwt — 21>, Jon—1 sin (20 — 1)wt.6 
n=1 n=1 
We multiply this expression by e*‘ and integrate term by term with 


5 Cf. G. N. Watson, Theory of Bessel Functions, (Cambridge, 1922), p. 22. 
6 In what follows, and here, also, the symbol J, shall indicate the value of the mth 


Bessel coefficient for the argument ib (i=4/—1). 
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respect to ¢ between the limits o and ¢. In this way we find that, as 
t is allowed to become infinite, /(¢) approaches 


a COS 2nwt + 2nw sin 2nwt 
a? + (21) ?w? 
asin (2 — 1)wt — (2m — 1)wcos (2n — a 


a? + (2n — 1)w? 





df oO 
I*(t) = pene) 4. as Fes 
a n=1 
(8) 





io 2) 
as 21>, Jen-1 
n=] 


The factor ein’ is now replaced by an expression similar to (7) in 
which 6 is replaced by —b. These two series are multiplied term by 
term and the corresponding harmonics are regrouped with the result 
that we find 


fo] 


J 2 re) re) JF iF on 
I*(t) = H > a ect + >| ( z. iss ) cos 2a 


I= 27 + 2m? = ny L\ 1. a? + (2 + 1)? w? 


2 2(2n + DwT J on+i ‘i 
+ (x aig ater =e, sin 2nwt 


— 210J J on+1-1 
" ( ) sin (2” — t 
» a? + (an +1 — 1)%w? ( I)w 


2 2i(2n +1 — 1)oJ Jonti1 ii 
_ t : 
+ ea Pile i ) cos (2n — 1)w 


That each Fourier coefficient is real is, of course, obvious. That it 
converges follows from the following considerations. We shall prove 
it for the coefficient of sin2mwt, neglecting the factor 2. Bearing in 
mind that | J,| =| J_1| , 


(2n + DoS J on+1 
a? + (2n + 1)*w? 








(9) 








ine} 





= 2, | Fhe 


l=—o 


2n 00 i) 
- bee ay +> | JF 120 | + > \dde | 
1=0 1=0 


l=2n+1 


( D\U/ b\ 22-2 b\)/ b\ 22 
< ss Qn (-) (-) = (—) = 
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1=0 l'\(2n — 1)! l=2n4+1 l\(2n -f 1)! 
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1=0 (an + 1)! } 





7 Cf. Watson, loc. cit., p. 16, (4). 
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b\2n 
< eb? /2(n+1) (2n +. 1) b\2” b \—2n sis (—) b2/4 | ° 
(n!)2 ( ) t ( ) me (2n)! ‘ 


It will follow, moreover, from this last expression that the Fourier 
series (g) converges. From the nature of the function J(#) and conse- 
quently also of J*(¢) its Fourier development is unique. It can be 
shown that J*(t) satisfies (1) and that equations (6) are satisfied by 
the corresponding coefficients of (9). In the course of this part of the 
work, use is made of the results in Watson (l.c.) on page 31, equations 
(3) and (4). 
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TRUE GROUND MOTION FROM MECHANICAL 
SEISMOGRAPH RECORDS 


A. B. BRYAN? 


ABSTRACT 


Records of several types of ground motion are obtained with two vertical mechani- 
cal seismographs of widely different characteristics, mounted on the same base and 
recorded on the same moving photographic paper. Companion records obtained are 
shown in one case to yield substantially the same calculated ground motion although 
the original records are very dissimilar. 

The record obtained with a seismograph is always affected to some 
extent by the characteristics of the instrument itself and never shows 
the true motion of the ground upon which the instrument rests. How- 
ever, if the instrument is a simple mechanical one of known character- 
istics, the true ground motion can readily be obtained from the orig- 
inal record by a simple but somewhat laborious procedure. For an 
electrical instrument of the type ordinarily used in geophysical ex- 
ploration the problem would be much more difficult. 

The problem of obtaining true ground motion has been considered 
in several published papers. H. A. Wilson? has discussed the problem 
thoroughly and has obtained the ground motions corresponding to 
some mechanical seismograph refraction records. More recently Karl 
Dyk? has described experiments in which records were obtained with 
a mechanical instrument mounted on a shaking table. The calculated 
motion of the table, as obtained from the record, agreed well with 
the directly observed motion of the table. In the present work, records 
are obtained simultaneously with two mechanical instruments of 
widely different characteristics mounted on the same base, and it is 
shown that although the two original records are not at all alike 
they yield the same calculated ground motion. 

The method of calculating true ground motion will be briefly re- 
viewed. The equation of motion of a simple mechanical seismograph 
is: 


(z) M(#+ 9) + R&+ Sx =0 
where x=displacement of the seismograph mass with respect to the 
instrument. 


1 The Carter Oil Company, Tulsa, Oklahoma. 
2H. A. Wilson, Physics, 12, 186, 1932. 
3 Karl Dyk, Bulletin of the Seismological Society of America, 25, 119, 1935. 
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y = displacement of the instrument and the ground on which it 
rests. 
M=mass of the suspended system. 
R=damping constant (viscous damping assumed). 
S=spring constant. 
From equation (1) it follows that: 


(2) —H= H+ (R/M)E + (S/M)x, 


(3) —jy=a+(R/M)x+ (s/M) [ xdt ++ Ki, 


t t t 
—y=n4t+ (R/M) { xdt + (s/m) f au f xdt + Kit + Ka, 
(4) 0 0 0 
Taking x=y=y=%X=0 at =o as initial conditions, then Ki= Ke=0 
and 


(Ss) -—y=ut (R/M) [dt + (s/a) fat fat. 


The instrumental constants (R/M) and (S/M) may readily be 
evaluated from a record taken with the instrument swinging freely. It 
may be shown that: 


(6) (R/M) = (2/T) loge r 
(7) (S/M) = (4m? + log.? r)(1/T)? 


where r=ratio of two consecutive swings in the same direction. 
T =period of oscillation of the instrument. 


In equation 5 the ground displacement (y) is expressed as the sum 
of three terms. The first term (x) is given directly as a function of 
(¢) by the original record, assuming that the mechanical and optical 
magnifications are known. The integral in the second term may be 
evaluated for any value of (¢) by graphically computing the area 
under the («) curve of the original record from its beginning up to 
the time (#). A second term may thus be obtained whose ordinate at 
any time is the area under the (x) curve from the beginning of the 
record up to that time. Similarly, the integral in the third term may 
be evaluated for any value of (¢) by graphically computing the area 
under the second term curve from the beginning up to that time, and 
a third term curve may be plotted whose ordinate at any time is 
equal to the area under the second term curve from the beginning up 
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to that time. The ground motion (y) is then obtained by adding to- 
gether the (x) curve of the original record and the second and third 
term curves, the latter two being multiplied by the previously de- 
termined instrumental constants (R/M) and (S/M) respectively. 
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In Fig. 1 are shown several pairs of records of the same ground 
motion made with two mechanical instruments of different character- 
istics. The two instruments were mounted on the same base and re- 
corded on the same strip of photographic paper. The No. 1 instru- 
ment had a period of 2.36 seconds while for the No. 2 instrument the 
period was only .455 second. Both instruments were very little 
damped, requiring fifteen or twenty seconds to come to rest when 
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disturbed. It will be noted that for ground motions obtained with 
explosives at short distances the records obtained in the two in- 
struments are very much alike although the natural period of the 
instrument does appear on the latter part of record 2B. For slower 
motions from other sources the two records are not at all alike. 
The original record 2D is shown as curve (x) in Fig. 2. The 
optical and mechanical magnifications of the instruments were not 
determined. The ordinates of all curves shown are of course multi- 
plied by this constant. Curve A, showing the area under the (x) curve, 
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was obtained by tracing the original record on coordinate paper and 
going through a process equivalent to counting squares under the 
curve. In this process the horizontal scale must be converted to time 
units, and variations in paper speed must be taken into account. It 
was convenient to use one-fiftieth second as the time unit. The period 
of the instrument was accordingly expressed in fiftieths of a second in 
all calculations. Curve A’ shows the area under the A curve, ob- 
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tained in a similar manner. It will be noted that curves A and A’ 
drift downward and become entirely negative at the end of the record. 
A drift of this sort is caused by a very slight error in placing the zero 
axis when the original record is traced or by a slight downward drift 
of the zero position of the instrument. 

The dotted curve in Fig. 3 shows the sum of the (x), A and A’ 
curves of Fig. 2, the latter two being multiplied by the values of 
R/M and S/M respectively, and should represent (—¥y) where (y) is 
the true ground motion. It will be observed that on account of the 
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previously mentioned drift the general trend of the curve at first 
rises and then falls rapidly to large negative values. Accordingly, a 
correction curve, shown as a Solid line, is arbitrarily drawn to follow 
this general trend as closely as possible and is subtracted from the 
dotted curve to obtain finally the true ground motion, as shown in the 
upper curve of Fig. 4. The other dotted curve in Fig. 4 shows the 
same ground motion as calculated in the same way from the com- 
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panion record 1D. In this case no correction curve was used. The con- 
stants of this instrument are R/M =.0116, S/M =.00287 and T=118. 
It will be noted that except for minor irregularities and differences in 
amplitude due to different magnification factors, the two instru- 
ments yield the same calculated ground motion. The original records 
are also for convenience reproduced in Fig. 4. It is seen that the 
No. 1 instrument in this case gives the true ground motion directly 
and that for this instrument no calculation is necessary. 

The principal difficulty with this method of calculating true ground 
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motion, aside from the labor involved, lies in the necessity for using 
an arbitrary correction such as that shown in Fig. 3. The assump- 
tion is made that very slow motions of the type shown by the correc- 
tion curve do not exist and hence should be subtracted as a correction. 
This assumption would for some types of motion be invalid. How- 
ever, if the motions are known to be reasonably rapid the correction 
does not introduce any uncertainty. A correction of the sort shown 
is made necessary by a Very slight error in drawing in the zero line on 
the original record. It may be shown that a correction curve ap- 
proximately the same as that shown in Fig. 3 would be made neces- 
sary if the zero line were drawn 0.1 mm too low at the beginning of the 
record and 0.6 mm too high at the end, and drafting errors of this 
magnitude are difficult to avoid. 

This work was done at the Rice Institute, Houston, Texas, with 
instruments loaned by the Humble Oil and Refining Company. The 
writer is obligated to Mr. N. N. Zirbel for assistance with the experi- 
ments and the computations. 














AN EXPERIMENTAL STUDY OF THE 
ELASTIC PROPERTIES OF ROCKS! 


J. M. IDE? 


ABSTRACT 


Laboratory determinations of Young’s modulus, rigidity, and compressibility were 
made on a set of representative rock samples. These measurements are compared with 
the theoretical relations between the elastic constants of an isotropic medium. Com- 
puted longitudinal and transverse wave velocities are compared with direct field 
measurements by Leet in granite and norite. Agreement is within 5 per cent for norite, 
and 20 per cent for granite. 


The interpretation of observed elastic wave velocities in the rocks 
of the earth’s crust is our principal basis for inferring its structure at 
inaccessible depths. The practical geophysicist is interested in time 
differences for waves travelling along various paths, rather than in 
the actual velocities in different rocks. However, the actual velocities 
are of fundamental importance when we deal with broader problems, 
such as the determination of the types of rock existing at the depths 
through which earthquake shocks travel. As Jeffreys remarks in his 
book, ‘‘The Earth”: ‘“‘The velocities of elastic waves in the upper 
layers of the crust give some information about the elastic properties 
of the materials. If this information is compared with experimental 
determinations of the properties of the rocks known from geology to 
be probable, we can proceed a long way towards identifying the main 
constituents of the crust.” 

As a basis for reasoning on geophysical problems, it is particularly 
important to know how well we are justified by experiment in apply- 
ing to large rock masses the usual assumptions underlying the theory 
of elasticity. 

If we make the customary assumption that a given elastic medium 
is homogeneous, isotropic, infinitely extended, and that Hooke’s law 
holds, then the following well known relations can be derived: 


(r) EB + 3E/u = 9 
i 
(2) ea: 2 Se 
ep 3—E/u 


1 The research described in this paper was performed at Cruft Laboratory, Harvard 
University, Cambridge, Massachusetts. 
? Shell Petroleum Corporation. 
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(3) v= — 


where £ is Young’s modulus of elasticity. 

B is the cubic compressibility of the medium. 

wis the modulus of rigidity (or shear). 

V is the velocity of a longitudinal elastic wave in the medium. 

v is the velocity of a transverse wave. 

p is the density of the medium. 
Several other relations can be written in terms of the above quantities 
and a, which is Poisson’s ratio of lateral contraction to longitudinal 
extension. Since these formulas are interdependent, the three equa- 
tions written above will be sufficient for the comparison of theory and 
experiment. 

In this research E and yu were experimentally determined for cylin- 
drical cores of rock, two inches in diameter and ten inches long. The 
compressibility, 8, had been previously measured by Dr. Zisman, of 
Harvard, on smaller cores of the same rocks. From these data we can 
test formula (1). 

For three rock formations, the Sudbury norite, Rockport granite, 
and Quincy granite, Dr. L. D. Leet of Harvard University*® has 
measured time-distance curves, using equipment similar to that em- 
ployed in seismic prospecting. From these curves experimental values 
of V and v were determined. Our measurements on rock cores from 
these formations enable us to compare Leet’s values with wave 
velocities computed from Young’s modulus and the modulus of 
rigidity, using formulas (2) and (3). 

The method of measuring Young’s modulus for the rock cores was 
devised by the author, and is described in detail in the Review of 
Scientific Instruments, October, 1935.4 The specimen, with metal foil 
cemented to its lower end, is placed upright on a thick steel disk. 
The foil and steel surfaces are separated by a thin sheet of mica. 
An alternating voltage of variable frequency applied across the elec- 
trical condenser thus formed, causes the rock to vibrate by electro- 
static traction. When the frequency of the voltage resonates with the 
natural longitudinal frequency of the specimen, the latter vibrates 
with increased amplitude, and radiates sound. This is detected by ear 
or by a piezo-electric crystal receiver. From the fundamental longi- 


* Physics 4, 375, (1933). 
4 See also J. M. Ide, Proc. Nat. Acad. Sci., February, 1936. 
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tudinal resonance frequency, fo, Young’s modulus is computed from 
the formula: 
1 /E 
| Sale” | ae 


al p 


where / is the length of the specimen, and p its density. The method 
is simple and rapid, and accurate within a few tenths of one per cent. 

Since the modulus of rigidity is related to the natural frequency of 
torsional vibration by a formula exactly analogous to this, it can be 
measured by an arrangement for setting the rock cylinders into tor- 
sional vibration. The method of this research is to fasten a small steel 
strip rigidly in a saw cut at one end of the rock sample, and produce 
torsion of the strip by two small electromagnets carrying alternating 
current from a variable frequency oscillator. The electromagnets are 
set on opposite sides of the strip, and produce a torque on it which is 
transmitted to the end of the rock. A similar device at the opposite 
end of the sample, with the electromagnets connected to an ampli- 
fier, is used to detect the vibrations of the cylinder. This arrangement 
is a modification of one used by Wegel and Walther of the Bell Labo- 
ratories,> in a study of the vibrations of metal rods. 

The compressibilities of these rocks were measured and the meas- 
urements described in detail by W. A. Zisman,® using a hydraulic 
testing machine to compress the specimen, while the yielding of 
the rock was measured by means of a sensitive rocking mirror piezo- 
meter. Petrographic descriptions of the rock samples may be found 
in Zisman’s papers. 

The results of our measurements, as applied to formula (1), are 
summarized in Table I. The table gives average values for three or 
more samples of each type of rock. The percentage deviation from 
the formula is given to the nearest per cent. 

The table shows that the theoretical formula which connects £, 
6 and up fits the experimental data with varying degrees of accuracy. 
For copper and steel, and for the diabase,—a compact fine-grained 
igneous rock, there is agreement between theory and experiment to 
about one per cent. This is excellent for measurements of elastic 
properties. For rocks which are less compact than the diabase, more 
porous, or composed of larger crystals, the formula does not hold so 


5 Wegel and Walther, Physics, 6, 141, (1935). 
6 Zisman, Proc. Nat. Acad. Sci. 19, 653 (1933); 19, 666 (1933). 
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well. The norite, quartzitic sandstone, marble, and granite show disa- 
greement increasing from three per cent to nineteen per cent. The 
Quincy granite is very coarse-grained, and the random variations in 
elastic properties from one sample to another are relatively large 
(fifteen per cent). 











TABLE I 
E 
EXPERIMENTAL TEST OF FoRMULA EB+3—=9 
Me 
E 
Substance EB+3— Deviation from o 
Me 
Steel 8.88 — I 
Copper 9-13 + 1% 
Vinal Haven Diabase Q.12 + 1% 
Sudbury Norite 8.76 — 3% 
Quartzitic Sandstone 8.15 — 9% 
Marble 7.48 —17% 
Rockport Granite 7.88 —12% 
Quincy Granite Feat —19% 








Another conclusion which followsfrom the data is that the disagree- 
ment between theory and experiment is always in the same direction, 
that is: the experimental value of E8+3E/y is uniformly smaller 
than we should expect. The quantity £@ is in general less than twenty 
per cent of the sum of the two terms of the formula, so that the 
quantity 3E/u is the one which shows the largest disagreement with 
the theory. 

We conclude then, that formula (1) agrees with experiment for 
metals and for a compact diabase rock, but that experimental values 
of E8+3E/u may be from nine to nineteen per cent lower than the 
formula predicts for such rocks as marble, quartzitic sandstone, and 
granite. We can also make a reasonable guess that disagreement be- 
tween theory and experiment will be still larger for the less compact 
sedimentary rocks such as shales, conglomerates and soft sandstones. 

In order to test formulas (2) and (3) by experiment, we compute 
longitudinal and transverse wave velocities from our measured values 
of elasticity and rigidity. We then compare these computed wave 
velocities with those derived by Leet’ from time-distance curves 
taken in the field. This comparison is made in Table II. 

It isevident from the table that for the Sudbury norite the formulas 


7 Op. cit. 
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agree with experiment to within about five per cent, while for the 
granites the figure must be raised to nineteen per cent. This result 
checks closely the comparison in Table I, where the agreement for 
norite and granite is three per cent and nineteen per cent respectively. 
Thus each of the three formulas from elastic theory is checked by 
experiment to about the same degree of accuracy. This accuracy is 
high for diabase and norite, and relatively low for granite. 


TABLE II 


COMPARISON OF FIELD AND LABORATORY DETERMINATIONS OF 
ELAstic WAVE VELOCITIES IN Rocks 











Transverse Longitudinal Differences 











Rock Wave Velocity | Wave Velocity 
v km./sec. v km./sec. Trans. Longit. 
— ag Kida | 8.22 Fila | — 2% | + 5% 
Rockport Granite a ae . - Su + 1% +19% 
a area | toon | ~18% | +14% 

















The reliability of the laboratory methods employed in this research 
is established by the excellent agreement with theory obtained for 
steel, copper, and diabase. Thus we must conclude that the dis- 
crepancies found for other materials are real, and that such rocks are 
therefore imperfect elastic media. The formulas of elastic theory are 
derived on the condition that the medium is homogeneous, and 
possesses the same properties in all directions. It is certainly reason- 
able that some types of rock should fit these conditions better than 
others. Our experiments show this to be true, and give the sign and 
order of magnitude of the discrepancies to be expected. 

It is possible that if the elastic constants of rock samples were 
measured at the pressures which exist even at a depth of a few hun- 
dred feet in the ground, theory and experiment might agree more 
closely. It is possible also that dispersion, or a variation of velocity 
with frequency, exists in rocks to a sufficient degree to explain part 
of the discrepancy between laboratory and field measurements. The 
laboratory experiments on elasticity and rigidity were carried out at 
frequencies near 10,000 cycles per second, while Dr. Leet measured 
the group velocities of transient shocks. 
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These questions are vital to the solution of important problems in 
geophysics, but they can be answered only by further experimental 
work on the elastic properties of rock, under the conditions which 
exist at great depths in the earth. 

In conclusion, the experiments described in this report suggest 
caution in applying the theory of elasticity to the surface rocks of 
the earth’s crust, whether they be sedimentary or granitic in char- 
acter. On the other hand, the results for diabase rock show that the 
theory can be expected to hold rigorously for the fine-grained igneous 
rocks which are presumed to lie beneath the oceans and the conti- 
nental granites. 














INTERSTITIAL WATER SATURATION IN THE 
PORE SPACE OF OIL RESERVOIRS! 


JAMES A. LEWIS? anp WILLIAM L. HORNER? 


ABSTRACT 


The presence of water in oil reservoir cores has been observed for several years by 
analysts using a direct method for fluid content determinations. The actual amount of 
water that is present in the sand is of importance in the estimation of oil reserves as 
well as in determining the type of fluid which a particular stratum will produce under 
conditions favorable for flow into the well. The core analysis profile of a well cored 
with oil (using Baker Core Barrel) is compared with three other wells which were cored 
with water. The analysis of the former core is in close agreement with respect to oil and 
water saturations, 30.58%, 44.40% respectively, with other cores in the area, and points 
to the conclusion that water was distributed throughout the sand section with no in- 
dication of gravitational separation at the time of coring. Coring with oil or the addi- 
tion of indicators to the fluid used during coring and the subsequent analysis of same 
are methods for determining correct interstitial water saturation. 


INTRODUCTION 


The presence of water in oil reservoir strata has been observed for 
the past several years by analysts using a direct method of analysis 
to determine the fluid saturation of oil well cores. The origin of this 
water saturation has been the subject of many discussions, and at 
present sufficient information is at hand to show that the data of 
per cent water saturation should be given weighty consideration 
when determining whether a well will produce water or commercial 
oil. It will be shown in the discussion which follows that water is 
present in the interstices of some oil sands at the time of coring and, 
depending on the actual amount of water present, is produced simul- 
taneously with the oil or is held in the pore spaces throughout and 
after the time of commercial operation. Water content data also are 
necessary in making quantitative estimates of the space available in 
pools for the accumulation and storage of oil. Calculations of fluid 
flow and estimates of per cent recovery and original oil content of 
reservoirs have long been made by assuming the complete oc- 
cupancy of the pores by hydrocarbons. This assumption seems logic- 
al, for it is often found that the oil strata yield pipe line oil without 
noticeable amounts of free water. However, the presence of water in 
the pores materially reduces the space available for the movement of 
and occupancy by oil, and consequently this paper may be of in- 


1 Read before the Society of Petroleum Geophysicists at Tulsa, March 20, 1936. 
2 Core Laboratories, Inc., Dallas, Texas 


353 








354 J. A. LEWIS AND W. L. HORNER 


terest to investigators of methods for predicting recovery and sub- 
surface fluid flow. It is interesting also to speculate as to the effect 
of this interstitial water saturation on the electrical methods in use 
at the present time in locating oil bearing horizons. The relative uni- 
form distribution of this water content throughout a sand of given 
equal characteristics, as pointed out later, as well as the fact that 
this water is in the position of contact in films of variable thickness 
with the sand, suggests to us the possibility of relatively low re- 
sistivities. In other words, it should be significant that the pores of 
many oil pools, instead of being filled entirely with oil, are filled 10 
to 50 per cent of their volumes with electrolytic solutions. Absolute 
resistivity values for oil sands in situ should therefore be lower than 
values for dry or artificially oil-saturated rocks. The presence of this 
water saturation in oil reservoirs is of importance in the evaluation of 
properties for secondary recovery methods as well as for those under 
natural production, and should be given consideration in the com- 
pletion of every well. The presence of interstitial water in an oil- 
bearing formation introduces the proposition that a definite oil-water 
contact does not exist in most reservoirs, and that water production 
in a well is not necessarily derived from a definite water strata or 
structural position, but can be produced from the same strata as is 
the oil. Geologic time has permitted all of the oil to rise above the 
water table in most homogenous sands, but has not let all of the 
water drop out of the oil zones. The contact usually referred to is, 
therefore, not a plane where one hundred per cent oil saturation over- 
lies one hundred per cent water saturation, but is a plane below which 
no oil production is obtained. 


ORIGIN AND SIGNIFICANCE OF WATER SATURATION 


The origin of this interstitial water in oil producing strata is open 
to question. However, it seems reasonable to us that, where the water 
percentage is not excessively high for sands of given characteristics, 
it is either connate water or water derived from the connate water 
which was present in reservoirs before the accumulation of oil, and 
because of the preferential wetting of silica by water in comparison 
to oil, is held in smaller capillaries and in openings where a condition 
of equilibrium exists between the quantities of water present and the 
diameter of the capillaries. As pointed out by Garrison (1) “In order 
that a well shall produce nothing but oil, it is necessary that the oil 
flow to the well bore through the larger pores at such a rate that it 
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does not cause the water to advance from its equilibrium position in 
the smaller capillaries to a point where it is forced along with the oil.” 

This condition of equilibrium represents varying percentages of 
water saturation in sands of different physical characteristics. Those 
sands which are formed of small grains will contain a higher percent- 
age of water in the state of equilibrium than will those with larger 
average grain size. Likewise, the displacement of water by oil will 
be more effective during the replacement period in the more perme- 
able formations, and will result in higher oil contents being found in 
these strata. The physical characteristics of a sandstone may be de- 
termined by the use of apparatus (2) (3) such as is found in petroleum 
production and research laboratories at the present time, and these 
results, when combined with determinations of fluid contents, yield 
data which indicate that there is a critical water saturation value for 
each and every sand that is worthy of note. 


DIRECT METHOD TO DETERMINE WATER SATURATION 


The presence of water in oil sand cores has been observed by the 
authors in each of several hundred determinations, and attention has 
been directed to investigating the question of its source both by 
direct and indirect methods. As is customary practice in the Ap- 
palachian Fields, practically all cores have been taken with a water 
head of from 60 to 80 feet in the hole. (Baker Method.) This practice 
has led to the general belief that the water which was found in all 
saturation samples was the result of penetration of coring fluid, and 
that water did not occur in a natural condition within the pores of 
the sand. Subsequent investigation and a change in procedure, 
whereby oil was used during the coring operation instead of water, 
facilitated the direct determination of water content, and at the same 
time presented an opportunity to compare the oil absorption and per 
cent oil saturation with previous results on cores taken with water. 

The first four core analysis profiles included herewith were drawn 
from laboratory results of cores from four wells located at the corners 
of a rectangle approximately 380 feet by 480 feet on the sides. C.W. 
2 and C.W. 3 as well as C.W. 4 and C.W. 5 are the pairs of wells 
which are spaced 380 feet with C.W. 2 and C.W. 5 on one 480 foot 
side, and C.W. 3 and C.W. 4 on the other 480 foot side. The order 
of coring was that designated by the numbers for the wells, with C.W. 
2, C.W. 3 and C.W. 5 being cored with water, and C.W. 4 being 
cored with oil. The area in which these cores were taken shows very 
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little structural relief on top of the sand, and within the distances 
represented by these wells, is of only general significance. A line across 
the diagonal of the rectangle from C.W. 4 to C.W. 2 follows a line 
from the structural low to the structural high. The inconsistency in 
total sand thickness is the result of an irregular base for the bottom 
of the sand. The natural reservoir pressure in the sand does not exceed 





Avg. Permeability, 


Depth, Feet Feet Sand Millidarcys 
45-8-56.2 10.0 63-4 
56.2-71.9 E57 28.8 
71-9-77-3 5-4 14.9 


Fic. 1. Core Analysis Profile C.W. 2. Cored with Water. 


15 lbs. per square inch gauge. The oil and water content of the sand 
was determined by distillation of fluids from the consolidated sample, 
and is the method which is in general use in the Appalachian Field 
for the analysis of oil field cores. By this method of analysis, the oil 
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and water content of a sample are determined simultaneously, and 
the results have been proven to such an extent that the proper evalua- 
tion of a property for expected recovery is a matter of routine an- 
alysis. 

Reference to Figs. 1, 2, 3, and 4 show permeability, porosity and 
oil saturation figures which are characteristic of this area. With the 





Avg. Permeability, 


Depth, Feet Feet Sand Millidarcys 
52-5-560.6 3-9 38.4 
56.6-71.3 14.7 35.0 
71.3-72-8 E<5 22.2 


Fic. 2. Core Analysis Profile C.W. 3. Cored with Water. 


exception of certain small thicknesses of sand which frequently are 
of a conglomeratic nature, the distribution of permeability values is 
such as to show a uniform sand with uniform flow characteristics. 
The results for these wells are summarized in Table I. The core in 
each instance has been divided into three sections on the basis of 
sand structure; the upper and lower divisions including in most cases 
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the conglomeratic strata which are considered too irregular for profit- 
able operation. The middle division, section #2, includes the sections 
of fine-grained sand which contain the higher percentage oil satura- 
tion values, and in turn, the highest oil contents. The bottom section 





Avg. Permeability, 


Depth, Feet Feet Sand Millidarcys 
31-0-36.5 4-9 77-9 
36.5-62.7 26.2 49.0 
62.7-68.1 inert 217.0 


Fic. 3. Core Analysis Profile C.W. 5. Cored with Water. 


of sand in C.W. 5, Fig. 3, possesses high permeability, with low oil 
saturation and high water saturation, and is a section which would 
act as a good path of travel for fluid entering or leaving the hole 
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should it be left open in the completed well. In this case, under equal 
pressure differentials, the bottom 4.5 feet of formation has a produc- 
tive capacity equal to the middle 26.2 feet of sand, and under natural 
production, under equal pressures, would produce a quantity of 





Avg. Permeability, 
Depth, Feet, Feet Sand Millidarcys 
26.0-34.6 8.6 16.8 
34.6-60.4 25e% 35-6 


Fic. 4. Core Analysis Profile C.W. 4. Cored With Oil. 


water equal to or greater than the amount of oil produced from the 
middle division. It is obvious that it is important to know the sec- 
tions of sand in a well that are producing a particular type of fluid 
in order that the proper completion of a well for the exclusion of water 
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may be known. It is likewise necessary and valuable to determine the 
net sand thickness which may produce oil under a given pressure at a 
definite rate, in order more thoroughly to evaluate the real worth of 
a productive zone. 

The core analysis profiles of the three wells just referred to were 
obtained by using water as the coring fluid. C.W. 4, Fig. 4, in 
the immediate vicinity of the previous wells was cored with oil from 
the top of the sand to a point 16.5 feet in, where water was used to 
core one screw. During the time that the well was cored with oil, the 
hole was bailed down after each three foot screw, and there was no 
evidence of water being present in the cavity. After taking the one 
screw with water in the hole, all fluids were then removed, and the 
coring continued with oil. The detailed laboratory results in Fig. 4 
represent permeability determinations at the points of the circles with 
fluid content determinations represented by solid lines for approxi- 
mately each foot of sand. All values are direct, and the actual values 
are shown on the upper line. 

The permeability values for this well are in agreement with the 
specific and average values of the previous wells and the distribution is 
normal. It is apparent that even though this well was cored with oil, 
the oil saturation of the core was in close agreement with the previous 
wells. Also it is apparent that the per cent water saturation is distrib- 
uted from top to bottom of the core in all sections of the oil zone and 
gives no evidence of gravitational separation. In addition to the above 
similarities between this core and the ones previously mentioned, 
another point is the per cent total fluid saturation in this sand, which 
is comparable to that present in the former cores. With a condition 
of 75% total fluid saturation with the water, oil and gas present in a 
state of equilibrium, and with an equal pressure around the sand 
biscuit at the time of coring, a condition suitable for fluid absorption 
where the Baker Method is being used is far beyond reasonable expec- 
tation. 

Inspection of the summary sheet gives the calculated average 
saturations for these four cores within the limits 28.07 to 30.65% oil; 
39.56 to 49.05% water, and 67.63 to 79.04% total fluid saturation. 
In particular, the well which was cored with oil had an oil saturation 
of 30.58% (not the upper limit of oil saturation values for these cores), 
and a water content of 44.40% saturation, which also is between the 
limits of maximum variation. 

Fig. 5 is the core analysis profile of a well where the water satura- 
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tion is higher than the amount found in any of the previous wells. 
With an average water saturation of 51.52%, and the low oil satura- 
tion of 21.97%, it is significant that water was the fluid that was 


& 





Avg. Permeability, 
Depth, Feet Feet Sand Millidarcys 
49-.I-5I.1 2.0 23.1 
i ar et 24.2 26.3 


Fic. 5. Core Analysis Profile C.W. 6. Cored with Water. 


bailed from this well. Under long continued operation, this well 
would not be a commercial producer, and its completion is not justi- 
fied. 
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OIL AND WATER SATURATION VS. PRODUCTION 


Sands of different physical characteristics have critical water 
saturation values that are controlled by the texture of the sand, and 
whenever the actual value exceeds this limit, water production can 
be expected; and whether or not oil is produced depends on the rela- 
tive amounts of both fluids that are present and their relationship 
to the total pore space. 

The results which follow are for cores in the Bradford, Pennsy]l- 
vania field on which production records are available. The exact 
water production in each case has not been recorded but is noted 
when the amount becomes greater than 3% of the total fluid pumped 
(Table IT). 


TABLE II 
FLUID SATURATION AND PRODUCTION FOR BRADFORD, PA., CoRES 
% Water % Oil Water 
Core No. Saturation Saturation Production* 
I 22 49 No 
2 28 39 No 
3 28 46 No 
4 29 48 No 
5 30 45 No 
6 33 56 No 
7 34 50 No 
8 34 4I No 
9 34 29 No 
10 35 35 No 
II 36 38 No 
12 36 30 No 
13 39 48 No 
14 39 40 No 
15 39 32 No 
16 40 45 Yes 
17 40 35 No 
18 40 24 Yes 
19 4! 38 Yes 
20 42 42 Yes 
21 43 42 No 
22 43 39 No 
23 43 17 Yes 
24 45 4I Yes 
25 46 48 Yes 
26 46 26 Yes 
27 Cr 38 Yes 
28 57 28 Yes 
29 59 27 Yes 


* Noted when water production exceeded 3% of the total fluid production during 
the first week on pump. 


INDIRECT METHOD OF WATER DETERMINATION 


Although oil may be used during coring to secure satisfactory 
samples for interstitial water determinations, the method has cer- 











364 J. A. LEWIS AND W. L. HORNER 


tain inherent difficulties which can be eliminated by the method of 
adding certain chemicals to the drilling fluid and the subsequent an- 
alysis of the water from the sand samples. This procedure (4) makes 
possible the division of the water content into that which was natural 
and that which was derived at the time of coring. Tests of this nature 
on Bradford Sand have shown that the amount of penetration is very 
small, and that for all practical purposes is of little value. As more 
permeable sands are cored which in their natural condition possess a 
lower total fluid saturation than does Bradford Sand, this question of 
contamination should be investigated more thoroughly. The character 
of the mud, the pressure differential which exists between the mud 
and the sand, as well as the contact area and the permeability, are 
factors which control the amount of penetration of the drilling fluid, 
and each should be given due consideration before a correct inter- 
pretation can be applied to interstitial water saturation. 


CONCLUSIONS 


From these data the following statements are given in concluding: 

1. Natural interstitial water of varying degree of saturation has 
been found to exist in certain oil reservoirs at the time of coring. The 
degree of saturation depends on the character of the sample as well 
as the distance the well is located from the advancing front of the 
edge-water encroachment. 

2. The interstitial water is distributed throughout the complete 
sand section, and gives no indication of gravitational separation. 

3. These data indicate that a definite oil-water contact does not 
exist in all reservoirs, and that water production in a well is not 
necessarily derived from a definite water strata or structural position, 
but can be produced from the same strata as is the oil. 

4. In addition to the usual factors in calculating productivity, the 
specific fluid saturation values, as well as the ratio of water and oil, 
are the controlling factors in determining the commercial value of a 


well. 
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SOME ASPECTS OF MULTIPLE RECORDING 
IN SEISMIC PROSPECTING 


PAUL W. KLIPSCH* 


ABSTRACT 

Considerable attention has been given to the use of more than one geophone on 
each recording channel with the hope of increasing the signal-to-noise ratio, where 
“signal” is taken to mean “‘recorded reflection’”’ and ‘‘noise’”? means any undesired re- 
corded amplitude. The average expectancy of this gain is evaluated analytically, and 
the statistical distribution examined. Since at different instants this gain possesses 
different values dependent on a large number of causes, these causes are assumed to be 
random in nature; a justification of this and other assumptions is explained, and the 
distribution of gain in terms of frequency of occurrence is expressed in the form of a 
probability curve. From this analysis an answer can be inferred to the question, “if ” 
is the number of geophones, what is the probability that a gain of k- +/x will occur?” 
The average expectancy of this gain is shown to be 


a 2 
G=—- V/n. 
VT 
This leads to the conclusion that a large number of geophones per channel is necessary 
to obtain consistently worthwhile values of gain in signal-to-noise ratio. 


Aspects incident to application are discussed, and sources of noise which are not 
random are mentioned. 


INTRODUCTION 


Multiple recording, or the use of more than one geophone per 
recording channel, has for its object an increase in the ratio of re- 
flected energy to other disturbances. Refracted energy, wind noise, 
seismic unrest, and uncontroled vibrations from traffic or other causes 
all tend to mask the ground motion which results from reflection 
from a deeply buried layer. 

The subject of multiple recording has been treated by Born! using 
Fourier analysis to show that certain components of the ground 
motion could be eliminated from the recording by using destructive 
interference, much in the manner of directive arrays of radio anten- 
nae. 

The recording may be regarded ascomposed of several components; 
the signal or reflection, regular noise, and irregular noise. The regular 
noise is that resulting directly from the shot; irregular noise includes 
seismic unrest and wind noise. The proposed method had for its 
purpose the elimination of certain components of the regular noise. 

There is considerable justification for an assumption that all noise 


* Independent Exploration Company. 
1 Born, W. T., ‘Notes on the Use of Multiple Shot Holes and Multiple Geophones,” 
presented at the Dallas convention, Soc. Petroleum Geophysicists, November 29, 1935. 
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be treated as irregular or random in its nature. If this be done, it is 
necessary to attack the problem from a different viewpoint involving 
the statistical approach.” It is the purpose of this paper to justify 
this basis for the study of the problem, and to carry out an analysis 
based on probability. 

In order to develop the idea that the ‘‘ground roll’ or refracted 
energy is random in its nature, it is necessary to review some of the 
fundamentals involved. This leads naturally to a study of the nature 
of the original disturbance, the paths through which the disturbance 
is propagated, and the manner in which superposition takes place. 


THE NATURE OF AN ARTIFICIAL SEISMIC DISTURBANCE 


The impulse generated by an explosion may be considered as a 
damped wave train. Considerable dispute arises as to the exact form 
of the wave sent out, but it should be safe to assume that the impulse 
gives rise to a ground motion which consists of a band or spectrum 
of frequencies, each (frequency) component of which consists of a 
damped wave train. Such a component can conveniently be repre- 
sented by a curve such as Fig. 1. Some components are probably 


: = nama 
TIME 
Fic. 1 








much more highly damped than others; it would be expected that all 
such components differ in amplitude, frequency, and damping, and 
suffer different degrees of attenuation in transmission. 

Fig. 1 can also be taken to represent an amplitude-space curve of 
a particular component of the traveling wave; in fact it is thé travel- 
ing wave in space which, passing a given point, gives rise to the ampli- 
tude-time curve appropriate to that point. 


SUPERPOSITION RESULTING FROM MULTIPLE PATHS 


The various ways by which elastic waves can travel from the shot 
point to the recorder provide an almost infinite number of motions 


* Klipsch, P. W., a prepared discussion following Mr. Born’s paper. See footnote 1. 
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at the recorder. In Fig. 2, a three layer earth is represented, in which 
the longitudinal velocities are indicated. In addition to the longi- 
tudinal velocites, there are the transverse and rotational velocities. 

Evidently, the impulse or wave train generated by the explosion 
will be propagated over various paths at a large number of velocities, 
arriving at the recording station as a series of superimposed wave 
trains, displaced in time, distorted in frequency characteristics, and 
altered as to magnitudes. The summation of these wave trains gives 
rise to a very complicated motion of the ground at any recorder, the 
complexity rendering an analytical attack practically impossible. A 
practical attack still remains, however, in a treatment of the com- 
plicated motion as random in nature. 

Si G. 
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Fic. 2 


Almost any event can be regarded as having a cause. For example, 
the flip of a coin shows heads. This is a result of its previous history, 
and the number of times it turned over. The causes are regarded as 
undetermined. In spite of the existence of causes, undeterminable in 
extent, the coin “‘obeys”’ the laws of probability. Similarly, in the case 
of a seismic refraction shot, the ground motion can not be determined 
a priori; the causes are undetermined until after the shot. Experience 
shows that no two records are alike, and even after the shot has been 
fired and the record analysed, it is impossible to learn all there is to 
know about the causal differences between any two locations. Enough 
independent causes, occurring at random, exist to justify the statistical 
approach. If an event, in the form of a particular ground motion, 
takes place at a point at a certain time, it will occur at a nearby point 
at some other time, the difference being expressible in terms of a 
phase angle of some particular frequency component. Even if the 
distance to the nearby point is not chosen at random, all phase angles 
within a limited range can be deemed equally likely. Several equi- 
distant colinear points will experience the event at different times, 
and, excepting first breaks or other regularly identifiable events, the 
phase angles of the event will be generally random in nature. All 
values of phase angles are thus equally likely. 
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One more important principle must be offered before the problem 
can be stated in a form capable of solution. This is that each fre- 
quency component, of which the disturbance is composed, can be 
treated separately, that all other frequency components obey similar 
laws, and that the total resultant is the summation of the separate 
components. 


GAIN IN SIGNAL-TO-NOISE RATIO 


In applying the principles of multiple recording, it is assumed that 
an array of shot holes or recorders can be sufficiently scattered so that 
the noise componenets add up out of phase, and at the same time 
that the array is sufficiently compact so that the signal, or reflection, 
adds up in phase. That this is possible depends on the relative ap- 
parent surface velocities of the signal and various components of the 
noise. 

Finding the signal-to-noise ratio gain thus depends on first finding 
how much noise will be reproduced by the array. If m is the number 
of geophones in the array, the signal amplitude will be ” times as 
great as for a single geophone. The noise, presumably approximately 
of equal amplitude at each geophone, but having different time or 
phase displacements, should produce a resultant less than m times 
the amplitude for a single geophone. The gain in signal-to-noise ratio 
would then be the ratio of gain in signal amplitude to gain in noise 
amplitude. The problem is then one of finding the noise amplitude as 
a function of x. 


ANALYSIS 


It was shown by Rayleigh? that the resultant of a large number of 
isoperiodic vibrations of arbitrary phase does not converge on any 
definite value. The average intensity is equal to m, where m is the 
number of sources of vibrations of unit amplitude, but for more than 
a superficial study it is necessary to inquire into the relative prob- 
abilities of resultants lying within certain assigned limits. It should 
be possible to state the probability that the resultant will be less 
than some specific value, say some function of x. 

The value of the probability arrived at by Rayleigh is 


(1) PH=r1—-erln 


3 Rayleigh, ‘On the Resultant of a Large Number of Vibrations of the Same Pitch 
and of Arbitrary Phase’’ Philos. Mag. X, 1880, page 73 (Scientific Papers, Vol. 1). 
See also “Theory of Sound” Vol. 1, pages 35-42. 
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where 
P=probability of a resultant amplitude less than r, 
r=resultant amplitude, r?=resultant intensity, 
n=number of unit vibrations of the same pitch and arbitrary 
phase. 
Equation (1) is plotted in Fig. 3, curve for n=, giving Pasa 
function of r?/n. 
AVERAGE VALUE OF RESULTANT 


The average value of the resultant may best be regarded as the 
expectation which is defined by Coolidge* 


(2) + pisi 
jul 


where ?; is the probability of occurrence of an amount s;. Thus the 
average value 


7= > rAP. 
Using the method of the Calculus, 


( 1 
3) A -{ r-adP 
0 
© 2r 
-{ e-?/n — dr 
0 nN 
- oe) ; rz 1/2 v2 
~ vi foerm(2)"{2) 
0 nN nN 


(4) = /a-—: 


It has already been shown’ that the average intensity 


(5) P= n. 
This may easily be verified by integrating the equation 
1 
(6) P= f r?-dP 
0 


4 Coolidge, Julian Lowell, ‘‘An Introduction to Mathematical Probability,’”’ Ox- 
ford, 1925. 
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ee 
= e—? /n — 2dr 
0 nN 


(7) =n. 
PROBABILITY OF A RESULTANT GREATER THAN AVERAGE 


To illustrate the practical implications of the statistical distribu- 
tion, consider the problem, “‘what is the probability of obtaining a 
resultant greater than the average?”’ The average resultant is given 
by (4) and we wish to find the probability that 7 exceed 7. The prob- 
ability, P, that the resultant is less than 7, is given graphically in 
Fig. 3 as a function of r?/n. Putting r=? gives (7)?/n=/4=0.785, 
for which value P=o.55, and 1—P=o0.45 which is the probability 
that the resultant is greater than the average value. 

Possibly of greater significance is the probability of a resultant 
greater than, say, one and one half times the average value. Here 
r=1.57, r?>/n=1.77, and 1-P=o0.16. Thus about 16% of the time 
the actual noise amplitude would be greater than 1.5 times the 
average noise amplitude. ° 


APPLICATION TO SIGNAL-TO-NOISE RATIO 


Assuming that several geophones on a single channel can be so 
placed that the signal, or reflection, arrives at all simultaneously, 
and that the noise arrives at random, it is possible to arrive at an 
expression for the expectancy of gain in signal-to-noise ratio. If &, 
denote the signal amplitude, X, the amplitude at each geophone, and 
n the number of geophones, then 


(8) é, — -/X,. 
The noise amplitude expectancy is a function of the noise arriving 
at each geophone and the number of geophones. If £; is the noise 


amplitude to be expected, or the average value of the noise amplitude, 
and X; the noise amplitude at each geophone, then, from (4), 


Vi 


 « 
ae ae 
(9) 2 
The signal-to-noise ratio with a single geophone per channel would be 
Xs 
(10) Ri = — 
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and the expected signal-to-noise ratio with geophones per channel 
would be 


(11) R, 


ms... 
E; 
The gain in signal-to-noise ratio resulting from changing from 1 to ” 
geophones per channel would have the expected value 
R, 





(12) : 
- avi. 


It should be remembered that the expected value of gain, G, is sub- 
ject to the same statistical variation as the average value of the result- 
ant, 7. That is to say, the probability of an actual gain, G, at some 
instant, less than G, is the same as the probability of a resultant, 7, 
greater than 7. Approximately half the time, the gain, G, will be less 
than G. 


9 
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SMALL VALUES OF 1” 


Rayleigh’s analysis (leading to the curve of Fig. 3) was based on 
the assumption that m is large. This immediately raises the question, 
‘chow large must m be so that the statistical analysis is at least ap- 
proximately correct?” 

The statistical distribution (similar to Fig. 3) has been worked out 
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for n=2 and n=3. Beyond this the problem becomes exceedingly 
difficult, involving finding the “volumes” of m dimensional spaces. 
Fortunately, however, a considerable convergence is indicated for 
the value of ~=3. These two cases will now be taken up in detail. 


DISTRIBUTION OF PROBABILITY OF RESULTANT OF TWO UNIT VECTORS 
Assume 2 unit vectors at angles @, and 42 with respect to an arbi- 
trary reference line. 
The resultant is 
r = [(cos 0; + cos 62)? + (sin 0; + sin 62)?]!/2 
(13) = [2+ 2 cos (0; — 62) ]!/2. 
By definition, all values of 6; and 42 are equally probable. Then the 
probability that the resultant is less than 7 is 
7 J frd0:d02 
14 ale’ sere 
I frd0,d0> 


where R represents a region of values of 6, and @2 within which the 
resultant is less than r and TJ represents a total region. Equation 
(14) is one of the fundamental equations of probability.® Since values 
of 0, and 62 greater than 27 represent repetitions it is sufficient to take 
T as the limits of 0; and #2 from o to 2m or from —7 to +7 in which 
case the denominator of (14) becomes simply 477’. 

The numerator of (14) is evaluated for different values of r by 


rearranging (13): 


2 cos (0; — 02) = r? — 2 


2 
(15) 6, = 02 + cos“! (= _ r). 


2 


The arc-cosine term has 2 values between o and 27 so the first integra- 
tion (on 6) is performed between the limits 


yz 
62 + cos7! (= _ :) 
2 


re 
0. + 27 — cos“! (= - r) 


2 


and 


5 Coolidge, loc. cit. page ro. 
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r 
62 — cos! (— —_ r). 
2 


The double integration in the numerator is equivalent to finding the 
area represented by the shaded part of Fig. 4a. To simplify the limits 
of integration Fig. 4a may be rearranged in the form of Fig. 4b in 


or simply 





LT 2r 


a 
9, 














° 
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) 
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Fic. 4 


which the areas are equivalent. Integration of 6. takes place between 


o and 27 so (14) becomes 
—1 (r?/2-1) 








2a 6,+¢ 8 
(16) P= So seen d0,d0. 
4m” 
2a r2 
Jo 2 cos“! ( — — x )dO. 
2 
4m” 
I r2 
T 2 


where P is the probability of obtaining a resultant less than r. Equa- 
tion (17) is plotted in Fig. 3, curve »=2. 


DISTRIBUTION OF PROBABILITY OF RESULTANT OF THREE UNIT VECTORS 
In the case of 3 unit vectors the resultant is 


r = [(cos 0; + cos 62 + cos 63)? 
+ (sin 0; + sin 02 + sin 63)?]!/2 
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(18) = [3 + 2 cos (0; — 62) + 2 cos (0; — 43) 
+ 2 cos (82 — 63) |!/2 


and the probability of a resultant less than r is 


‘“ Sf frd0,d0ed05 
Sf rd0\d02.d03 


The denominator of (19) is simply 82°. Since the numerator is not 
amenable to the analytic methods used in (14), resort was had to 
graphical methods. 

Equation (18), when studied geometrically, leads to a nest of 
cylinders; holding any one angle constant, as a parameter, gives a 
trace of the cylinders in a plane, and then permitting the parameter 
to vary generates a solid. This leads to the conclusion, which could 
have been arrived at intuitively, that one of the angles can be set 
equal to zero and the volume of a solid reduced to the area of a plane 
surface; this amounts to letting the arbitrary reference line follow 
one of the vectors. An analogy for the two vector cases is that in Fig. 
4b, the length of a line instead of the area of a surface could have 
been used as a measure of the probability. This amounts to saying 
that (13) reduces to 


(13b) r = (2+ 2 cos 6)!/? 





(19) 


and (14) reduces to 


p _ Ladd 


27 





(14b) 


where the R region is simply a length of line. In like manner, (18) 
reduces to 


(18b) r = [3 + 2 cos (01 — 62) + 2 cos 0 + 2 cos O2]!/? 


and (19) to 


(0b) tia j Sred01d82 


4m? 


Holding 7 constant in (18b) gives rise to a closed curve, the area of 
which, divided by 47’, is the probability that the resultant lies be- 
tween o and 7. 

Equation (18b) is plotted in Fig. 5. Measurement of areas results 
in the probability curve for n=3 in Fig. 3. 
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ANALYSIS FOR LARGE VALUES OF ” HOLDS FOR SMALL VALUES 


Since as ” increases, convergence must occur on the curve for large 
values of 2, and since the curve for n=3 is a fair fit to the curve for 
n=, it is to be inferred that for m>3 the curve for large values of n 
applies to a very close degree of approximation; that is, the statistical 
analysis for large values of ” applies closely for all values of ” greater 
than three. 


k 





-7r (6) 0; Tr 


FIc. 5 


PRACTICE 


Equation (12) shows that the average gain in signal-to-noise ratio 

is approximately 
G = 1.13Vn. 

Evidently, to obtain an average gain of between 11 and 12 would 
require 100 geophones per channel, or 600 for a 6 channel recorder. 
This assumes that other factors, such as depth to which shots and 
geophones are buried, are not altered. Overlooking the obvious ab- 
surdity, consider the economics of placing a large number of geo- 
phones, or shots, below the irregularly attenuating and retarding top- 
soil, in order that all factors other than the number of units may be 
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the same. Twelve hundred feet of geophone holes is obviously out of 
the question. 

On the other hand, suppose the cost be avoided by placing the 
geophones on the surface. Then the top-soil characteristics will fre- 
quently produce errors in reflection time and step-out as great as 
would be caused by ten times as many feet of ordinary weathering. 
The low-pass filter characteristics of a soft topsoil will also reduce or 
even eliminate the irregularities known as phase breaks by which a 
particular half cycle of a reflection is best identified for correlation 
purposes. Such loss in detail may more than offset the theoretical 
advantage. 

Equation (12) may be regarded as establishing a maximum gain to 
be expected from the compounding. An assumption on which the an- 
alysis is based was that a configuration is possible by which the signal 
is recorded in phase and the noise out of phase. It should be realized 
that such a configuration may not be possible in all cases. For ex- 
ample, a reflection from a dipping bed, shot down-dip, may have a 
very low apparent surface velocity, and an array spread which is suffi- 
ciently great to cause the noise to be random may receive the signal 
with such phase lags that the signal amplitude is less than m times 
that for a single geophone. In such a case, the gain would be materially 
less than given by (12). Incidentally, records shot up-dip should gen- 
erally be expected to display better gain in signal-to-noise ratio than 
records shot down-dip. 


PREVIOUS REFLECTIONS AS A SOURCE OF NOISE 


A type of noise so far not considered is that produced by several 
successive reflections occurring at such close intervals that they tend 
to mask each other. Suppose that a multiple shot or geophone system 
were so perfect that all noise of the irregular variety were completely 
eliminated. Whatever unwanted ground motion exists, of the nature 
of ground roll or refraction energy, is assumed to be reduced to a 
level below perceptibility. Only reflections would be recorded. If 
several reflections occur in quick succession, each reflection leaving 
behind it a decaying wave train, interference between reflections will 
occur which will tend to mask the information obtainable from each 
taken separately. That is, the decaying wave train of a reflection will 
constitute an interfering signal, or to all intents and purposes a noise, 
as far as the succeeding reflection is concerned. At least in the Gulf 
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Coast, where reflections appear on the records in great numbers, this 
would prove serious. 

It is to be inferred, then, that what is regarded asa signal at one 
instant becomes at some later time a noise. Being part of a signal, 
it cannot be eliminated by a system which merely discriminates 
against random noise. Evidently, the multiple system cannot produce 
any advantage in the solution of a problem of this kind. 
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Prospeccao Geophysica em Sdéo Paulo. By Irnack Carvalho do Amaral and Henrique 
Capper Alves de Sousa. Bulletin 10 of the Department of Mineral Production, 
Republic of Brazil, Rio de Janeiro, 1936. (In Portuguese.) 


Part I deals with the magnetic instruments, methods and anomalies found in the 
areas of Sao Pedro-Xarqueada. Askania vertical and horizontal intensity magnetom- 
eters, and stationary instruments for photographic registration of diurnal and other 
changes, were used. Most of the anomalies were found to be related to diabase sills, 
lacoliths and dikes. Depth calculations by the methods of Tiberg, Haanel and Thalen 
did not correspond closely with depths found in wells, but were of the same order of 
magnitude. Depth to basement rocks could not be calculated from the magnetic re- 
sults, and the authors suggest the use of seismic methods. Tables showing the magnetic 
properties of the rocks in the area are included. 

Part II deals with the theory, instruments, field procedure and results of gravita- 
tional surveys in the Sao Pedro area. Small Bamberg balances were used. Porosity and 
density of rocks in the area are listed. Calculations for interpretation are made on the 
assumption that masses have an infinite extent in one direction. Maps showing gra- 
dients, curvature values, and isogams are presented for the Sao Pedro and adjacent 
areas. 

Part III deals with a magnetic reconnaissance in the vicinity of Pitanga and Pau 
D’Alho, where anomalies are related to diabase sills, dikes, and lacoliths, and to struc- 
tural elevation of rocks of relatively high permeability. 

This publication of 102 pages contains little new in the way of theory. It will be of 
interest to those studying the geology of Brazil, and to those engaged in magnetic sur- 
veys involving intrusive and extrusive igneous rocks. Requests for the publication 
should be addressed to Bibliotheca do Departamento Nacional da Produccao Mineral, 
Avenida Pasteur 404, Praia Vernelha, Rio de Janeiro, Brasil. 

Joun H. WILson 


“Deformation of Rocks under High Confining Pressures. I. Experiments at Room 
Temperature.” Journal of Geology, Vol. xliv, no. 5, pp. 541-577. By David T. 
Griggs. (1936). 

This article describes experiments on the flow and rupture of rocks under compres- 
sion, tension, and torsion, while at the same time subjected to a high confining pressure 
supplied through a liquid surrounding the specimen. The hydrostatic pressure of this 
liquid could be measured very accurately and could be maintained constant. In addi- 
tion, a “differential” stress was applied to the specimen, and the deformation was 
measured directly. By using the high pressure technique of P. W. Bridgman the confin- 
ing pressure was carried up to 13,000 atmospheres, equivalent to a depth in the earth’s 
crust of 28 miles, and four times that available to F. W. Adams in his pioneering experi- 
ments (1901-1917). 

The experiments of Adams had demonstrated beautifully that under high confining 
pressures rocks behave plastically. However, he was unable to measure his pressures 
accurately, to hold them constant, or to rupture the rocks under a confining pressure 
and examine the fracture after the experiment. All these defects have been remedied in 
the present apparatus. 
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The results of the present experiments contradict the ancient geological dogma that 
under a sufficient confining pressure a rock would flow indefinitely like a viscous fluid. 
Instead, it was found that whatever the confining pressure, up to the limits available, 
a differential pressure great enough to cause flow, if maintained long enough, would 
always eventually produce rupture. Nor did rapid loading under these circumstances 
produce rupture at a less stress. On the contrary, the specimen ruptured at less stress, 
but with greater deformation, under slow loading. This is really to be expected from the 
first finding. 

The ultimate strength, and the amount of deformation preceding rupture, were 
vastly increased by a high confining pressure. But the elastic limit, and the behavior 
below that point, were little affected. Quartz remained brittle up to the highest pres- 
sures available. 

The article is illustrated by many beautiful photographs of ruptured or deformed 
specimens, and by photomicrographs of the rocks before and after deformation. 

It is to be hoped that the author’s plan to continue these experiments at higher 
temperatures, and at still higher pressures, will soon be carried out. 

M. Mott-SMItH 


“Shearing Phenomena at High Pressures of Possible Importance for Geology.” Journal 
of Geology, Vol. XLIV, pp. 653-669, 1936. By P. W. Bridgman. 


Contributions from Bridgman’s high-pressure laboratory at Harvard are always 
of interest. In the present article he discusses flow, and the shearing stresses necessary 
to produce flow under high pressure for a wide variety of substances. Rigidity, or the 
resistance to shearing stress, is one of the fundamental constants of the seismologist, 
and Bridgman here deals with shearing stresses under conditions of great pressure com- 
parable to those under which seismic data are obtained. The experiments were made at 
room temperatures, but at pressures as high as 50,000 kg./cm.,? corresponding to those 
prevailing at depths of 166 km. Few rocks now available for examination at the surface 
of the earth are likely to have been buried so deep. 

“‘The phenomena found proved to be both varied and unexpected.” Great pressure 
produced a many-fold increase in shearing strength, 10-fold for most metals, goo-fold 
(an extreme case) from paraffine. Even under conditions of extreme distortion crystal 
structure persisted and at high pressures polymorphic transitions occurred in many 
cases, some with notable volume changes. Repeated deformation increased shearing 
strength in some cases, decreased it in others. “‘In the case of graphite ‘‘all the lubricat- 
ing qualities associated with its initially easy internal slip are eventually lost, and at the 
end graphite becomes an abrasive, actually rubbing metal off the hardened steel sur- 
faces.” Plastic flow in forty per cent of the cases was accompanied by deformation, by 
internal rupture, and rehealing. Permanent, non-reversible changes in density occurred, 
and even chemical reactions. A number of substances showed important time effects. 

“The existence of these various time effects suggests crucial questions with regard 
to behavior in geologic times. We are here, of course, in a field of almost pure specula- 
tion; but there are examples of such abrupt freezing of internal mobilities with decreas- 
ing temperature or increasing pressure that I think we should not be frightened by the 
bogy of even geologic time, but in our speculations freely use such concepts as shearing 
strength, admitting that the apparent strength in 10,000,000 years may be a little lower 
than for laboratory experiments, or the strength for 100,000,000 years a little lower than 
for 10,000,000, and also recognizing that the difference between geological and labora- 








380 REVIEWS OF RECENT PUBLICATIONS 


tory behavior may be a strong function of the material. We must at any rate recognize 
that the existence of these time effects means that no numerical datum or descriptive 
concept formulated in simple language is sharply applicable to geological phenomena.” 
“Geology is rapidly approaching a point where an ultimate problem which has been 
staring it in the face can no longer be sidestepped, namely, to determine the actual 
specific physical and chemical behavior of those materials which do actually constitute 
the earth’s crust.” 

The reader must keep in mind the fact that all these experiments were conducted 
on very thin films of material, only a few thousandths of an inch thick, that the pres- 
sure applied was unilateral, and not the conditions of cubic compression and enhanced 
temperatures that exist in the depths of the earth. 

C. MAyNaArD Boos 


“Magnetic Investigations in Southwest Alabama” by J. Brian Eby and E. G. Nicar; 
Geological Survey of Alabama, Bulletin No. 43, University, Ala. (June, 1936). 


The report covers the results of field work with the Askania-Schmidt vertical 
magnetometer done by the junior author. The area worked includes all that part of 
Alabama lying within sixty-five miles of the Alabama-Mississippi boundary between 
the Gulf of Mexico and the northern portions of Choctaw, Clarke, and Monroe coun- 
ties. A map of that area, on a scale of three miles to the inch, showing anomaly values in 
gammas and locations of stations occupied, and contoured on an interval of twenty 
gammas, accompanies the report. 

A few clear concise paragraphs on the topography and surface geology of the area, 
are followed by a discussion of the principles involved in magnetometer work and their 
application to this particular area. By way of illustrating these principles, reference is 
made to igneous plugs in Arkansas and Mississippi, and some salt domes and anticlines, 
all previously mapped with a magnetometer by the junior author. It is stated that the 
magnetic relief of the igneous plugs was always positive, ranging from 600 to 2600 
gammas, that of salt domes negative, and that of other anticlines either positive or 
negative, the magnitude of relief in the last two classes not exceeding 44 gammas and 
often requiring correction for the effect of deeply buried igneous masses. The effect of 
paramagnetic Tertiary sediments in producing the magnetic pictures of these last two 
types is discussed. 

Six magnetic anomalies are pointed out as a result of the magnetic work in Ala- 
bama. Two of these are of the high-relief, positive type. It is stated that they can be 
interpreted as deep-seated plugs of igneous rock. While it is said that such plugs may be 
either erosional remnants on the surface of the basement or actual intrusions into the 
sediments, the authors neglect to point out that the plugs might also be basic intrusions 
or inclusions in an acid basement, which would not necessarily rise above the general 
level of the basement and thus might have no effect on the structure of the overlying 
sediments. 

The other four anomalies explicitly referred to are marked by closures of twenty to 
forty gammas or by a mere spreading of the contours. Several other such anomalies 
shown on the map are not mentioned in the text. The authors point out that on the 
flanks of the major magnetic highs one could scarcely expect to find minor closures 
without correcting for the rapid lateral increase of magnetic intensity due to the 
major feature. Finally, they state that the general magnetic picture apparently is not 
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affected by two previously known geological structures: the Hatchetigbee anticline and 
the Jackson fault. Eight well logs from the region are appended. 

Since the average density of magnetic stations read is only about one for each five 
square miles, the twenty gamma contour interval used appears excessively small to the 
reviewer. It has the effect of showing up small closures at intersections of the traverses 
run (on which stations were placed at one-mile intervals) and tends to give the im- 
pression of a featureless field within the polygons between the traverses, although 
doubtless several such small closures also exist within those polygons. The multiplicity 
of contours may somewhat obscure a rather interesting regional relationship shown by 
the map, namely, that the main streams of the region and their estuary in Mobile Bay 
follow in a general way the magnetic lows and distinctly avoid the major magnetic 
highs. This could be interpreted as meaning that the regional magnetic picture reflects 
in a rough way the configuration of the surface of the basement, and that the central 
drainage system of Alabama is geologically at least as old as early Cretaceous, occupy- 
ing a geosyncline with, in general, a north-south axis. 

JosepH L. ADLER 


“Some Magnetometric and Gravimetric Surveys in the Transvaal.’”’ By O. Weiss, 
D. T. Simpson, and G. L. Paver, Department of Mines of the Union of South 
Africa, Geological Series Bull. No. 7, 1936, pp. 27, plates 8. 


This short but interesting report covers test application of the magnetic and gravi- 
tational] methods of prospecting to geologic problems in seven areas in South Africa: 
1) Magnetic traverses across the Lower Witwatersrand System. Anomalies rang- 
ing between + and — 10,000 gammas, close relation between anomalies and ge- 

ology were observed. 

2) Torsion balance survey across the Doornkop Fault; a text book example of 
gradient and differential curvature anomaly. 

3) Torsion balance survey across the Witpoortje Fault and Lower Witwatersrand 
beds; an example of anomalies which are easy to interpret if you know the 
geology. 

4) Successful determination of the hade of the Witpoortje Fault by a magnetic 
survey. 

5) Torsion balance survey for locating nickel ore pipes. Brilliant indication is 
given of the pipes in spite of a regional gradient of more than 50E. 

6) A Magnetic survey over Gossan, indicating the usability of magnetic methods 
in prospecting for nickel ore pipes. 

7) Magnetic survey with negative results over chrome seams. 

This is a valuable collection of cases of the application of these two methods. 

D. C. BARTON 
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James A. Moore 

W. Harlan Taylor, T. L. Allen, O. S. Petty 
Elay Kenton Myers 

R. D. Wyckoff, Paul Weaver, E. A. Eckhardt 
Leo James Peters 

T. A. Manhart, Paul Weaver, E. A. Eckhardt 
Pedro Rey 

E. E. Rosaire, Derry H. Gardner, Louis Statham 
Terence Charles Richards 

J. H. Jones, John H. Wilson, A. L. Smith 
John Frederick Riddle 

J. S. Watt, R. R. Thompson, Louis Statham 
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Frederick Ernst Romberg 

L. G. Howell, W. W. Rust, Jr., J. S. Watt 
William Schriever 

J. S. Watt, R. R. Thompson, M. A. Arthur 
Marcel Schlumberger 

E. G. Leonarden, E. E. Rosaire, F. M. Kannenstine 
Ellis William Shuler, Jr. 

Eugene McDermott, J. C. Karcher, W. E. McDermott 
Frank Benjamin Smith 

O. S. Petty, T. L. Allen, W. Harlan Taylor 
John Barron Souther 

W. W. Rust, Jr., J. S. Watt, Louis Statham 
Melvin Carhart Terry 

J. S. Watt, W. W. Rust, Jr., W. B. Lewis 
William Edmund Wallis 

G. F. Kaufman, Dave P. Carlton, J. E. LaRue 
Waldo William Waring 

J. S. Watt, Leland L. Antes, A. Frosch 
Loren Elvin Whitehead 

L. F. Athy, L. A. Scholl, Jr., E. V. McCollum 
Joe Tracy Williams 

J. S. Watt, Leland L. Antes, A. Frosch 
Floyd Willis, Jr. 

M. C. Born, W. T. Born, J. E. Owen 
Robert Ellerston Woods 

Roy F. Bennett, A. L. Smith, John H. Wilson 
Hughes Mead Zenor 

J. S. Watt, W. W. Rust, Jr., W. D. Mounce 


FOR ASSOCIATE MEMBERSHIP 


Willis Oliver Bazhaw 

T. L. Allen, O. S. Petty, W. H. Taylor 
Richard Everett Blight 

E. P. Evans, O. C. Lester, L. Y. Faust 
Richard Bunting Coit 


Eugene McDermott, J. C. Karcher, J. E. Jonsson 
Jack F. Fore 

E. E. Rosaire, J. W. Flude, T. I. Harkins 
Lawrence Frank Guseman 

C. C. Zimmerman, L. A. Scholl, Jr., C. T. McAllister 
Audio Gray Harvey 

Eugene McDermott, J. E. Jonsson, Barney Fisher 
Hubert Slocum Harvey 

Derry H. Gardner, L. Atmar M. Barnette, J. S. Watt 
Jesse Duff Hatch 

Derry H. Gardner, L. Atmar M. Barnette, J. S. Watt 
Homer Lafayette Heggy 

W. R. Ransone, Barney Fisher, J. E. Jonsson 
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Walter Johnson 

Robert Saibara, E. E. Rosaire, J. W. Flude 
Cecil Klemme Jordan 

L. F. Athy, L. A. Scholl, Jr., E. W. McCollum 
Ernest August Kiesler 

Barney Fisher, Eugene McDermott, J. C. Karcher 
Kenneth Morrison Lawrence 

O. C. Lester, Jr., John L. Ferguson, L. Y. Faust 
Manuel Calva Lecuona 

J. M. Adler, H. W. Falkenhagen, J. M. Golden 
William Isaac Mayfield 

T. L. Allen, O. S. Petty, W. H. Taylor 
Marion Jefferson Moore 

T. L. Allen, O. S. Petty, W. H. Taylor 
Pedro Eugenio Narvarte 

. T. L. Allen, O. S. Petty, W. H. Taylor 

Robert George Nisle 

H. B. Peacock, Eugene McDermott, J. C. Karcher 
James Angelo Smith 

W. M. Rust, Jr., W. B. Lewis, J. S. Watt 
Elden Edward Stroshine 

W. Harlan Taylor, T. L. Allen, O. S. Petty 
Brunislaus Stephen Ulak 

John W. Flude, E. Stiles, C. Maynard Boos 
Daniel Rey Vercesi 

W. G. Green, T. A. Manhart, F. B. Leedy 
Joe LaFayette Whittenberg 

J. S. Watt, Leland L. Antes, A. Frosch 
Robert Edward Wilder, Jr. 

E. E. Rosaire, T. I. Harkins, J. L. Lake, Jr. 


EXECUTIVE ORDER NO. 21 
CHANGE OF STATUS 


By authority of Article III-B-6 of the constitution, and in accordance with their 
Executive Order No. 21, the Executive Committee has advanced the following from 
associate active to membership: 

B. O. Winkler 
D. F. Broussard 
EXECUTIVE ORDER NO. 18 
FOREIGN MEMBERSHIP FEES AND SUBSCRIPTIONS 


The Executive Committee has authorized the addition of fifty cents per year to the 
membership fee for foreign members and fifty cents per year to the subscription price 
to foreign addresses to take care of additional postage, loss in exchange, and handling 
charges incident to such memberships and subscriptions. Individual numbers will be 
billed to foreign addresses at regular prices plus twenty cents per number. 








